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SECTION 1 


PROGRAM OVERVIEW AND SUMMARY 


INTRODUCTION 

NASA is conducting a flight control research program in digital fly -by-wire technology 
using a modified F-8C aircraft. The first phase of this program used Apollo hardware 
to demonstrate the practicality of digital fly-by-wire in an actual test vehicle. For 
the second phase, conventional aircraft sensors and a large floating point digital 
computer are being utilized to test advanced control laws and redundancy concepts. 

For control law studies, the new computer's speed, memory, and floating point capability 
removes most of the constraints imposed by more typical flight control machines. 

As part of this research activity, Honeywell initiated work in 1974 under Contract 
NAS1-13383 to provide a system of digital adaptive flight control laws for flight test 
in Phase II. These control laws were to be adaptive in the sense of altering the control 
law on the basis of normal sensed information. One adaptive design was to be refined 
for flight test from a comparison of several candidate concepts. As part of the ground 
rules, measurements were restricted to rate gyros, accelerometers, and servo position. 
Air-data was excluded because aircraft like the F-8C, whose performance requirements 
are readily satisfied with air-data -scheduled control laws, benefit most directly from 
adaptive control through the elimination of air-data sensors. This will be particularly 
valuable for future aircraft requiring flight critical redundant control systems and 
sensors. The control laws were further constrained to be compatible with the existing 
airframe without structural modification. Hence, they use only existing elevator, 
rudder, and ailerons as control effectors— each powered by existing actuators. 

Three candidate control algorithms were successfully designed from distinct adaptive 
concepts: 

• Gain adjustment based on self-excited limit cycles, 

• Gain adjustment based on explicit identification with Liapunov-stable model 
trackers, and 

• Gain adjustment based on explicit identification with maximum likelihood 
estimation (MLE). 


channel implementation. Several Kalman filter channels operate at fixed locations in 
parameter space. Likelihood functions are computed for each. Sensitivity equations 
are then solved only for the maximum likelihood channel and used to interpolate from 
there to the correct parameter value. 

Five parallel channels suffice- to handle the F-8C aircraft over its entire operational 

flight envelope. They estimate three parameters: surface effectiveness (Mg e ), 

pitching moment due to angle-of-attack (M^), and airspeed (V). Estimation accuracy 

depends strongly on the signal levels in the control loop. For the small test signals 

tolerable in operational situations, errors of 10 to 20 percent in M and 20 to 30 percent 

in and V are typical in six-degree-of-freedom simulation runs. Theoretical accuracy 

analyses confirm these error levels. The gain adjustment in the pitch and lateral 

control laws is a function of estimated only. 

6 e 

CONCEPT COMPARISON 

The preceding three designs were compared on the basis of performance, growth 
potential, and computer requirements. An overview of the comparison is given in 
Table 1. The MLE design was selected primarily on the basis of growth potential. 


TABLE 1. OVERALL COMPARISON OF CONCEPTS 


Characteristic 

Model tracker 

High-gain 
limit cycle 

Maximum 

likelihood 

estimation 

Performance 

Acceptable to good 

Good 

Good to excellent 

Growth 

potential 

Low— limited to 
single variable 
explicit gain 
schedule 

Low — limited to 
single variable 
implicit 
adaptation 

High — multiple 
parameter gain 
adjustment possible 

Computer 

requirements 8 

<5 

Minimal 

Low 

Medium 


a 

These requirements are relative to A P-101 capacity. 


The entire control system is expected to fit comfortably in the F-8C's AP-101 flight 
computer. The control law plus adaptive algorithm (excluding redundancy management) 
are estimated to consume less than one-half of the frame time available on that marh<m» 
(20 msec for a 50 Hz update rate) and to require approximately 2500 memory locations. 

4 




SECTION 2 


INTRODUCTION 


BACKGROUND 

NASA is presently conducting research in digital fly-by-wire technology using a modified 
F-8C aircraft as the test vehicle. In 1972-73 during the first phase of this program, 
the F-8C was successfully flown using an Apollo guidance computer, an inertial measure- 
ment unit, an analog (electrical) backup system, and an electrohydraulic actuation 
system. The purpose of these flights was to demonstrate the practicality of a digital 
fly-by-wire flight control system in an actual test vehicle. 

The second phase of the research program supports technology development of advanced, 
reliable control systems for future aircraft. Primary objectives of the second phase 
are to demonstrate benefits of advanced control' laws and to study redundant. flight control 
mechanizations. The test aircraft will use a triplex AP- 101 digital computer and triply - 
redundant aircraft gyros and accelerometers. Dual air-data measurements of angle-of- 
attack, Mach number, and pressure altitude will be available as well as a single side- 
slip sensor. Of special importance is the extensive computing power provided by the 
Large, floating point computer in the test aircraft. This capability makes it possible . 
to flight Test advanced control laws without the constraints imposed by less capable 
flight control computers. 

To date, two sets of control laws have been developed for flight experimentation in 
Phase II: a "CCV Package" and an "Adaptive Package. " The first package is described 
in Reference 2. It provides basic command augmentation functions, outer loops, and 
control modes for various Control Configured Vehicle concepts applicable to fighter 
aircraft. The second control law package is described in this report. It consists of 
the same command augmentation functions found in the CCV package but adds explicit 
on-line identification capability for adaptive control. In addition to these packages, 
research has also proceeded on several other candidate adaptive algorithms for the _ 
F-8C 3 ’ 4 ’ 5 and on sensor fault detection and isolation algorithms. ’ Some of these 
may reach flight experimentation later in the program. 
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MOTIVATION FOR ADAPTIVE CONTROL 

NASA's primary reason for studying adaptive control laws for the F-8 aircraft is to 
identify potential performance and system effectiveness benefits offered by modern 
adaptation/identification concepts and to demonstrate these in an actual digital flight 
control environment. Relative to these objectives, adaptive control is viewed in this 
report as one of several competing ways to achieve successful control actions in the 
face of aircraft plant uncertainties. Whether it offers benefits for a particular design 
application then depends on two primary factors: 1) the nature and magnitude of plant 
uncertainties, smd 2) alternate methods available to achieve successful control. The 
first factor determines feasibility. Most adaptive techniques work well for some types 
of uncertainties but not for others. The second determines viability. Adaptive controls 
may or may not be competitive with alternate approaches on economic or operational 
grounds. 

For inner-loop stability and command augmentation of fighters, the feasibility factor 
tends to be positive. Plant uncertainties generally take the form of unknown parameters 
in an otherwise known model structure (i.e,, coefficients of linearized equations of 
motion). Their range of uncertainty is largely due to widely varying flight environments 
(dynamic pressure, Mach, angle-of-attack) and large configuration variations (c.g, 
j locations, fuel and payload, geometry). However, individual coefficients are strongly 
interrelated, and only a few must be known accurately in order to design good control 
laws. As a net result, the important uncertainties reduce to a small number of slowly 
varying parameters. They usually include M g , sometimes also M ff , and occasionally 
other surface effectiveness coefficients. Adaptive controls based either on explicit 
identification or on implicit control gain adjustment stand a good chance of success 
with these kinds of uncertainties. This is demonstrated by three separate F-8C adaptive 
designs andsimulator evaluations presented in this report and will be verified, hopefully, 
in upcoming F-8C adaptive flight tests. 

Two alternate approaches for dealing with uncertain aircraft dynamics are to use low- 
gain flight control laws which can tolerate the full range of uncertainties directly or 
to provide external air-data for gain adjustment. The first approach is generally ruled 
out on the basis of performance. Low-rgain systems capable of flying the entire flight 
envelope are unable to satisfy performance requirements for most of today's fighters 
and will be even less acceptable for future conditionally stable relaxed-static-stability 
aircraft. The second approach is more successful. . It exploits the strong correlations 

6 



which exist between externally measurable air-data variables (q, Mach, «) and unknown 
aircraft parameters. Thanks to these correlations, control law gains which are originally 
computed as functions of M fi , M^, etc, can be approximated in terms of q, Mach, and/or 
«. Although not as good as the original functions, such approximated "air-data schedules" 
are adequate to meet performance-requirements in most current flight control prpblems. 

Given these alternatives, the competition among flight control diesign techniques reduces 
to a tradeoff between the economics, reliability, and operational features of air-data- 
scheduled systems versus adaptive systems. Historically, this tradeoff has favored 
adaptive techniques because sufficiently accurate and reliable air-data sensors were 
difficult to build and maintain. In fact, the desire to avoid such sensors was the primary 
motivation for most adaptive control development efforts of the 1960's, including NASA’s 
X-15 research program, 6 * 8 * 9 * 10 * Ut 12 * 13 * 14 the F-101 flight tests and X-20 simula- 
tions, 16 the production F- 111, * 7 and several more recent flight tests with the 
F-4, Qf CO urse, the measurement issue also motivated sensor development 

efforts. These have produced improved air-data systems which have now reversed the 
tradeoff, making air-data scheduling the preferred approach in most current single- 
channel analog autopilots. 

The present status of the adaptive versus air-data trade-off is not expected to be permanent, 
however. There are three important trends which will affect future choices. The first 
of these is digital control. Flight control systems are now being mechanized in on-board 
digital 'Computers 2 * whose steadily decreasing price tags per function promise to reduce 
the costs of adaptation substantially. 

The second trend is redundancy. Future control configured fly-by-wire aircraft may 
need dual.^ triple, or even quad redundant flight control systems to insure safety of 
flight. Air-data-scheduled approaches will then suffer the added cost and complexity 
of multiple air-data sensing systems, and they will also face the difficult problem of 
finding suitable probe locations for adequate data quality and channel tracking (to permit 
effective redundancy management) but with enough dispersion to avoid common hazards. 


* All of these aircraft can be flown quite well with approximated air-data schedules, 
provided that accurate measurements are available. 
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The final trend is increased control performance. Several anticipated flight control 
applications appear to require control gain functions which cannot be adequately approxi- 
mated by air-data schedules. Prominent examples are active control of flexure 22 ’ 23 

24 

and coupled control modes. Successful flexure control laws, for instance, have 
been found to depend critically on accurate vehicle data: mode frequencies, mode shapes, 
rigid body/flexure/control interactions. Such data at structural frequencies are difficult 
to acquire. The same is true for coupled control modes which actuate several surfaces 
simultaneously in order to achieve deliberate coupling or decoupling of responses 
(e. g. , flight path changes without attitude changes). Here the relative surface effective- 
ness values and several other stability derivatives must be known with high accuracy. 

These flight control technology trends provide strong motivation for continued research 
and development in adaptive control. Studies with the Phase II F-8C are particularly 
appropriate because its redundant, digital, fly-by-wire configuration is a forerunner of 
the evolving technology. Of course, all concepts and design processes studied should 
be general enough to apply to other aircraft as well. This is important because the F-8C 
itself is not a state-of-the-art fighter nor is it difficult to control with standard air-data 
scheduling methods. 

CANDIDATE ADAPTIVE CONTROL LAWS FOR THE F-8C 

This report presents design details and simulator evaluation results for three candidate 
adaptive! control laws for the Phase n aircraft. Each design provides basic command 
augmentation functions using only inner-loop inertial sensors (gyros and accelerometers) 
and no air-data measurements. Hence, each is a potential competitor with air-data- 
scheduled methods. The purpose of these designs is to establish feasibility of adaptive 
controls in a digital fly-by-wire setting and to examine the relative merits of different 
adaptive techniques. They are not intended to demonstrate control performance superiority 
over air-data scheduling (the F-8C offers little hope here) or to show explicit economic, 
reliability, or operational advantages. The latter must wait for detailed hardware 
designs and trade-offs. The three designs achieve adaptation by combining existing 
control laws from the M CCV Package” with different on-line adaptive gain adjustment 
processes. The gain adjustment processes are: 

!♦ Explicit parameter identification via Liapunov-designed model tracking with 
gains stored as functions of parameters. 
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Explicit parameter identification via maximum likelihood estimation (MLE) 
with gains stored as functions of parameters, and 

3 # Direct gain adjustment to maintain is elf- excited limit cycles. 

These specific adaptive techniques were selected from the brief survey of applicable 
adaptive concepts summarized in Table 2, This survey categorizes various proposed 
adaptive methods into four groups : 

1, Optimal Adaptive 

These are concepts based on the "dual control" problem originally posed by 

OC * . . ' ( • 

Feldbaum. They involve optimal identification and control done simultaneously 
according to an overall criterion of goodness. This problem is so difficult 
to solve that to date only approximate solutions are available even for trivial 
cases, 

2. Performance Adaptive (implicit) 

These are concepts based on on-line measurement and adjustment of performance 
variables. Included are high-gain, model-following schemes and model reference 
systems. 

3, Separate Identification and Control (explicit) i r/<o 

These concepts are based on assumed separability of the identification function 
and the control function in adaptive control. While this has been shown to be 
subop tiirial in the "dual control 1 ' sense, the assumption seems reasonable in 
many cases* Table 2 lists several proposed approaches for the separated 
functions of identifying and recomputing the control law. 

4. Learning / Self- Or ganiz ing 

This group consists of concepts with "memory, " Once a control law has been 
found to work well for a recognizable environment, the law is stored in memory * 
and called for again when the same environment reappears. The "control law" 
in this case includes not only gain parameters, but the entire structure of the 
relationship between measurements and control commands. 
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TABLE 2. ADAPTIVE CONCEPTS 


[ Optimal adaptive 

Dual control approximations 

[26-29] 



Performance adaptive (implicit) 

High-gain model following 




► Limit cycle control 

[8-14] 



► Energy balance control 

[15,16] 



— ► Damping ratio control 

[17] 



► Test signal detection 

[18] 



Model reference 




Performance gradient design 

[30] 



Liapunov design 

[31-34] 



Separate identification and control (explicit) 

Identification 

■* 

Control 


Batch processors 


Algebraic methods 


Equation error least squares 

[35,36] 

Deadbeat control 

[35, 44] 

Output error least squares 

[37-40] 

Pole -zero placement 

[52] 

Maximum likelihood 
Instrumental variables 

[41-43] 

[44] 

Equation coefficient placement 
Optimization methods 

[32] 

Recursive processors ' 


Ricatti equation solutions 

[46,50] 

Least squares 

[45, 46] 

Single stage optimization 

[53, 54] 

Extended Kalman filtering 
Stochastic approximation 

[47] 
[48, 49] 

— ► Stored solutions 


■ » Model tracking 

[50,34,51,19,20] 

— ► Constant bandwidth 

[19,20,55] 

Learning /SMf- organizing 

Trainable controllers 

[56, 57] 



Learning systems 

[58] 





ggygPQ,! references sire provided under ench category to describe pur ti c ula r techniques 
more fully. These were selected to treat flight control applications where possible. 
However, they should be viewed as representative references only. Given the vast 
adaptive literature, no attempt was made to be all-inclusive. 

The surveyed concepts range from ideas which exist only on paper to systems already 
tested in flight. The latter are highlighted by solid arrows in the table. They include 
four high-gain, model-following designs: 

1 . 

2 . 

3. 

4. 

and two explicit model-tracking identification schemes with stored gain functions: 

19 

1. The "Navy Adaptable System" (flight tested on F-4 ), and 

20 . 

2. A simplified "SIDAC" system (flight tested on an F-4 ). 

Selections of candidate approaches from Table 2 were made to span a range of performance 
and complexity, to maximize flight-test value, and yet to bound the risk of achieving at 
least one flight-test-quality design. No concepts were selected from the Optimal Adaptive 
and Learning /Self-Organizing categories, for example, because these have only limited 
development status with respect to practical applications. A high-gain limit-cycle 
system was selected because it offers very low risk (versions have flown) yet still has 
flight-test value if the design can demonstrate that improved limit-cycle detection and 
tracking logic can overcome weaknesses of earlier analog systems. The maximum 
likelihood procedure was chosen as representative of recent modern estimation theoretic 
approaches which promise very good identification performance at the expense of on- 
board computer complexity, while the model -tracking approach was chosen to examine 
identification capability with more modest computer requirements. 


Limit-cycle control (operational on the X-15 ), 

15 

Energy balance control (flight tested on the F-101 and demonstrated for 
the X-20 on iron-bird simulations with triply-redundant control hardware ), 

17 

Damping ratio control (operational on the F-lll ), 

18 

Test signal detection (flight tested on an F-4 ) # 


* This is discussed in more detail in Section 10. 
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Successful simulator-quality designs were developed for these three selected adaptive 
concepts, and each was evaluated on NASA Langley's F-8C simulator. 59 As expected, 
control performance for all three was roughly equivalent. Identification capability 
turned out to be best for the MLE concept. Hence, this design was refined further and 
streamlined to fit into the AP-101 flight computer. It is now recommended for flight- 
test implementation. 

DOCUMENT ORGANIZATION 

This report is subdivided into 12 sections. The next section. Section 3, presents a 
list of symbols used throughout. Operational requirements and ground rules for the 
adaptive designs are set forth in Section 4; aircraft models and simulations are summarized 
in Section 5; and Section 6 presents basic theoretical identifiability limitations for the 
F-8C which are imposed by the ground rules of Section 4. The common CCV control 
law structure for all three concepts is then described in Section 7. Sections 8, 9 
and 10 are devoted to the individual design details for the maximum likelihood, model- 
tracking, and high-gain limit-cycle concepts, respectively. Concept comparisons and 
selection of the most promising concept for flight test are described in Section 11, and 
Section 12 presents the conclusions and recommendations. 
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SECTION 3 


SYMBOLS 

Operators 


Arg {minf(x)} 

X 

E(* ) 

minimizing argument of function f(«) 
mathematical expectation 

exp(* ) 

exponential 

ln (‘ ) 

natural logarithm 

tn det(* ) 

- natural log of determinant 

p(x/y) 

- conditional probability distribution of x given y 

S. 

Laplace operator 

z 

delay operator 

• 

11 

time derivative 

A(0 

increment 

v(*) 

gradient vector with respect to parameter vector £ 

v 2 (* ) 

second partial derivative matrix with respect to 
parameters £ 

V* 

p-th component of V(* ) 

s 

summation 

M 

absolute value 

H'llJ. 

T 

quadratic form x Mx 
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SYMBOLS (Continued) 


Superscripts 

A 

( ) 

n 


< > 


( >° 


(i) 


Subscripts 

u 


trm 


( ) 


m 


( )• 


WL 


<>, 


Upper Case Symbols 

A 

B 

By 

C 

c* 


estimated value 
one-step predicted value 
value for parallel channel i 
nominal value 

value at trimmed flight 

(1) measured value 

(2) model value 

water line value 
value at time t. 


discrete system dynamics matrix 

(1) discrete system input matrix 

(2) residual covariance matrix 

sensor bias 

discrete system constant input matrix 
response variable N z + V cq q 
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SYMBOLS (Continued) 


D 

(1) measurement matrix, y due to u 

(2) design parameter matrix in model tracker 
equations 

F 

- continuous system dynamics matrix 

V S) ’ V s >’ G W g (s) 

gust filter transfer functions 

G c *’ g lat' g roll’ g rss 

control gains 

G r g 2 

continuous system input matrices for u, T| 

H ‘ H x* H u* H W 

- measurement matrices 

I 

identity matrix 

J 

partial likelihood function L - 1/2 H -tndetB 

K. K* 

- continuous and discrete Kalman filter gains 

K p* K I’ k l* k 2 

- control gains 

K c 

- critical gain 

L 

- likelihood function 

L p » L r » L 6a* L 6r 

- rolling moment coefficients due to indicated variables 

L , L » L 
u v w 

- gust field scale lengths 

M 

- number of parallel channels 

M q » M a . M 6 (M 6e ). M. a 

- pitching moment coefficients due to indicated 

variables 

M 

o 

- trim pitching moment 

M 6t 

- true value of 
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SYMBOLS (Continued) 


(1) number of data samples 

(2) discrete sensor noise intensity matrix 

Yawing moment coefficients due to indicated 
variables 

normal and lateral acceleration 
weighting matrix in Liapunov function 
a priori parameter covariance matrix 
partitions of P 
probability 

continuous sensor noise intensity matrix 
transformation matrix 

thresholds in model tracker error nonlinearity 
sequence of N control inputs 

(1) velocity 

(2) Liapunov function 

true airspeed 

crossover velocity in C* response 

lateral force coefficients due to indicated variables 

sequence of N measurements 

normal force coefficients due to indicated variables 



SYMBOLS (Continued) 


Lower Case Symbols 



hysteresis parameter 

parameter vector with components c^ 

a priori estimate of £ 

true value of £ 

partitioned parameter vector 

true value of £^» £2 

sensor displacement from c.g. (4.62m) 
model tracker error vector 
nonlinear model tracker error 
ij-th element of matrix F 
gravity 

disturbance process to represent gravity variations 
in maneuvers 

ij-th element of matrix G 

altitude 

index of the minimum-L channel 
roll rate 

dynamic pressure 

(1) yaw rate 

(2) number of measurements 
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t 

u 


u 


e 


U t 

u g 


w 


w 


Greek Symbols 
Upper Case 


r 

$ 

o 


SYMBOLS (Continued) 
time 

(1) Control input vector 

(2) forward velocity perturbation 

elevator servo command 
test input 

forward gust component 

(1) lateral velocity perturbation 

(2) error parameter in model tracker equations 

lateral gust component 
vertical velocity perturbation 
vertical gust component 

(1) forward position perturbation 

(2) state vector 

(1) lateral position perturbation 

(2) measurement vector 

parameter errors in model tracker equations 


discrete system noise input matrix 
tracker error dynamics matrix 
subset of parameter space 
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VL, VL 2 



SYMBOLS (Continued) 

- partitioned gradient vector 

partitioned second partial derivative matrix 


Greek Lower Case 


g 


6 ( 6 ), 6 , 6 , 6 ,, 6 . 

e or r rt 


e, «j, «2 


V Ae k 


D 

e 




xy 


angle - of - att a c k 

gust angle -of- attack 

angle of sideslip 

gust angle of sideslip 

aerodynamic surface positions 

error quantities 

generic likelihood filter states 

( 1 ) damping ratio 

(2) dummy argument for values of parameter 
vectors 

white noise process 
pitch attitude 

parameter adjustment gains in model tracker 

(1) exp(-At/ t) 

(2) design parameter in model tracker 
Kalman filter residuals 

white noise process 
correlation coefficient 



SYMBOLS (Concluded) 


parameter in Liapunov function 
standard deviation of variable x 

(1) time constant 

(2) data length 

roll attitude 
yaw attitude 
natural frequency 

design parameters in model tracker 



SECTION 4 

DESIGN GROUND RULES AND REQUIREMENTS 


GROUND RULES 

The aim of the design effort reported here was not to develop new theoretical procedures 
or algorithms but to turn existing concepts into flight-worthy control laws for the specific 
test aircraft. The concepts and design processes, however, should be general enough 
to apply to other aircraft as well. Other more specific ground rules imposed on the 
design are described in the following paragraphs. 

Inputs 

The adaptive system must operate in the presence of normal pilot inputs and also when 
such inpats are absent. Any test signals required for the latter case must be small 
enough not to interfere with the aircraft's mission. This generally means that test 
input normal accelerations, as sensed at the pilot station, should be below (0.2 - 
! 0.3 m/sec 2 ), and lateral acceleration should be even lower in the range of (0.1 m/sec ) 
This ground rule establishes a crucial distinction between identifiers designed for 
operational adaptive controls and those designed for post flight-test data processing 
applications. 42, 43 In the latter case, test inputs are deliberately large and often 
optimized for identification accuracy. 60, 61 For operational adaptive controls, the 
emphasis must be to make these inputs small and, hopefully, unnoticeable. 

| Sensors 

| The system was constrained to operate with the aircraft's existing surface (servo) 

1 position transducers and inertial sensors (rate gyros and accelerometers), but without 
attitude measurements and without air-data measurements (dynamic pressure, velocity, 
angle-of-attack). These limitations are imposed by the sensor complement available 
on the test aircraft and, in the case of air data, by the motivation in Section 2. 



Computer Capacity and Sample Rate 


Since control law calculations typically consume only a small fraction of the total 
computational load in a flight computer (the majority is I/O, self test, mode, and 
redundancy control), the adaptive law's goals were restricted to 25 percent of the available 
frame time of the AP-101 and to a "reasonable" allocation of memory. The sample 
rate of the computer was pre-specified to be 50 Hz. This gives five msec real-time 
per sample in which to complete all control and identification calculations. The sample 
rate is also high enough (for the F-8C) to produce no substantial differences between 
direct digital design (discrete-time control laws designed for discrete-time models) 
and continuous time design with after-the-fact discretization. 

Control Surfaces and Actuators 

The standard elevator, aileron, and rudder surfaces were available for adaptive control. 
Leading edge and trailing edge flaps were assumed to remain undeflected or to follow their 
open-loop scheduled minimum drag positions. Actuation systems for these surfaces 
were represented by models incorporating projected characteristics of the improved 
systems being Installed on the test aircraft. 

OPERATIONAL REQUIREMENTS 

While no detailed specifications were written for the adaptive control laws, all algorithms 
were required to operate satisfactorily under the following conditions: 

1. Periods of quiescence 

2. Periods of high turbulence 

3. Large (all-attitude) maneuvers 

/ * 

4. Landings and takeoffs 

Satisfactory operation was demonstrated by simulator evaluations under these conditions. 


Satisfactory operation for these mission phases was assured by inhibiting the adaptive 
process and using fixed gains. 
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PERFORMANCE CRITERIA 


The perfor m ance criteria deliberately concentrated on identifier performance* This 
was done since overall closed-loop performance for the F-8C provides a less sensitive 
measure of adaptive capability.'- Four "measures of goodness" were used: 1) identifica- 
tion accuracy in steady flight, 2) convergence characteristics, 3) tracking characteristics 
for standard flight transition, and 4) responses to major maneuvers and configuration 
changes. 

Accuracy at Fixed Flight Condition 

Identification accuracy (and thus control gain accuracy) were verified by subjecting 
the closed-loop adaptive system to a standard sequence of test conditions while in 
trimmed flight at fixed flight conditions which span the F-8C envelope. The test sequence 
consisted of a period of quiet, followed by a 6 m/sec 2 (20 ft /sec 2 ) pilot doublet C* -command, 
followed by a period of atmospheric turbulence, and then again by a 6 m/sec 2 doublet 
This entire sequence was repeated with and without sensor noise. The standard 
sequence was used to simplify data taking and comparison between concepts. 

Convergence Characteristics 

Convergence with identified parameters initially in error by ± 12 db was checked during 
both gust disturbances and pilot commands. Acceptable performance is required within 
three seconds after initialization, with the disturbance and system engagement occurring 
simultaneously. 

Varying Flight Conditions 

Time-varying flight profiles were used to exercise adaptation capabilities. These 
consisted of acceleration and deceleration maneuvers without pilot stick activity. Since 
all concepts emphasized adaptation to Mg g , the profiles were chosen to encounter large 

variation in Mg . 

6 

Major Maneuvers 

The closed-loop adaptive systems were also tested with large maneuvers, including 
rolling maneuvers, high-g turns, and step disturbances in angle-of-attack and slideslip. 
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SECTION 5 

MODELING AND SIMULATIONS 

AIRCRAFT MODELS 

The aircraft models used for the design work were derived from nonlinear aerodynamic 

data for the F-8C supplied by NASA, Linear airplane models were obtained by numerical 

■ 2 I 

differentiation techniques as described in the CCV control law final report. The 
resulting pitch and lateral linear models have the form 

x = Fx + G.^u + G 2 w (1) 

y = H x + H u + H w 
J x u w 

The structure and numerical values of these matrices are given in Reference 2, The 
state, x, is defined in Table 3. The control input, u, is given in Table 4, and the 
components of the disturbance, w, and output vector, y, are listed in Tables 5 and 6, 
respectively. Twenty-five flight conditions used for design are listed in Table 7, 

The aircraft modeling was completed by adding actuator dynamics for the control 
surfaces, including the gust model, and by defining sensor noise. 


TABLE 3. STATE DEFINITION 


Symbol 

Definition 


P 

Roll rate (rad/sec) 


r 

Yaw rate (rad/sec) 


V 

Lateral velocity (m/sec) | 

Lateral 

9 

Bank angle (rad) 

axis ; 

* 

Yaw angle (rad) 


y 

Lateral displacement (m) 


q 

Pitch rate (rad/sec) 


w 

Vertical velocity (m/sec) 


u 

Forward velocity (m/ sec) 

Longitudinal 

9 

Pitch angle (rad) 

axis 

h 

Vertical displacement (m) 


X 

Forward displacement (m) 
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TABLE 4* CONTROL DEFINITION 


Symbol 

.Definition 


6 a . 1 . 

Differential aileron deflection (rad) 

Lateral 

*r 

Rudder deflection (rad) 

axis 

®rt 

Differential elevator deflection (rad) 


6 e 

Symmetric elevator deflection (rad) 


®sb 

Speed brake (rad) 

Longitudinal 

8 g 

Gear (%) 

axis 

*t .. 

Throttle <%) 



TABLE 5. DISTURBANCE VECTOR 


Symbol 

Definition 


u g 

Forward gust (m/sec) 


v g 

Side gust (m/sec) 


W g 

Vertical gust ( m/sec) 


TABLE 6. OUTPUT VECTOR 


Symbol 

Definition 

p 

Roll rate (rad/sec) 


q 

Pitch rate (rad/sec) 


r 

Yaw rate (rad/sec) 


N z 

2 

Normal acceleration (m/sec ) 

N 

Lateral acceleration(m/sec^) 

1 


_j 
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TABLE 7. TWENTY -FIVE FLIGHT CONDITIONS 


1 

- 1 h 

Mach 

q(psf) 

' 

q(n/m 2 ) 

V (m/sec) 

“trm (de S> 

Condition 

(Kft) 

(Km) 

1 

20 


6.1 

0.670 

305 

14,603 

212 

3.45 

Cruise 

20 


6,1 

0.670 

305 

14,603 

212 

6.10 

‘ 

AN = lg (climb) 
z 

20 


6.1 

0.670 

305 

14,603 

212 

12.12 

AN = 3g (climb) 
z 

20 


6.1 

0.670 

305 

14,603 

212 

4.32 

Cruise 

20 


6.1 

0.400 

109 

5,219 

126 

8.86 

Cruise 

20 ; 

6.1 

0.900 

551 

26, 382 

285 

2. 18 

Cruise 

40 


12.2 

0.700 

134 

6, 416 

207 

6.73 

Cruise 

i 

40, 

12.2 

1.200 

395 

18,913 

354 

2.72 

Cruise 

io; 

3 

0.800 

652 

31, 218 

263 

1.96 

Cruise 

0 


0 

0.700 

725 

34, 713 

238 

1. 86 

Cruise 

0 


0 

0.300 

133 

6, 368 

102 

7.64 

Cruise 

0 


0 

0.530 

416 

19, 918 

180 

2.88 

Cruise 

i ' 

20| 

6.1 

0.600 

245 

11, 731 

190 

4.25 

Cruise 

20 


6.1 

0. 800 

435 

20, 828 

253 

2.54 

Cruise 

40 


12.2 

0. 800 

175 

8, 379 

236 

5.15 

Cruise 

4or 

12.2 

0.900 

222 

10, 629 

266 

4.08 

Cruise 

0 


0 

0.189 

53 

2,538 

64 

7.48 

Power approach 

0 


0 

0.219 

71 

3,400 

75 

2.76 

Power approach 

20; 

6.1 

0. 670 

305 

14, 603 

212 

2.12 

AN = 0. 5g (dive) 
z 

20! 

6.1 

0.600 

245 

11, 731 

190 

15.45 

AN * 3g (climb) 
z 

40 

12.2 

1.400 

537 

25, 712 

414 

2.64 ' 

Cruise 

40 


12.2 

1.300 

463 

22, 168 

384 

2. 65 

Cruise 

30 


9.1 

1.200 

633 

30,308 

364 

1. 92 

Cruise 

30 

9.1 

1. 100 

532 | 

25,472 

334 

1. 89 

Cruise 

30' 

9.1 

0.600 

158 

7,565 

i 

182 

6.09 

Cruise 

















ACTUATOR MODELS 


Actuator models consisting of a first-order lag for the primary actuator and a second-order 
lag for the secondary actuator were included for the elevator aileron and rudder. The para- 
meters were again taken from Reference 2 and are listed in Table 8, The secondary actuator 
dynamics were approximated as a unity transfer function for the design of the model tracker 
and the maximum likelihood estimation concept. The hysteresis values shown were estimates 
of the values expected in the Phase II hardware following modifications . The complete model 
is shown in Figure 2 . 

For the majority of the simulation runs, the hysteresis was set to zero. This was considered 
since the adaptive algorithms used measured secondary servo position, and in flight the 
hysteresis of the primary actuator is expected to be small. This results from the aerodynamic 
loading of the surface. 


TABLE 8. ACTUATOR MODEL PARAMETERS 


Actuator 

Primary (first-order lag) 

Secondary (second-order lag) 

Pitch 

= 0.0800 
2a = 0.00436 rad 
Rate limit = 0.44 rad/sec 

a) = 125. 6 rad/sec 
£ = 0.7 

2a = 0. 00175 rad 

Roll 

t 2 = 0. 0333 
2a =0.006 

Rate limit = 2.44 rad /sec 

Same as above 

Yaw 

Tg = 0.0400 

2a = 0.00855 rad 

Rate limit =122 rad/sec 

Same as above 


GUST MODEL 

The gust model used has the form attributed to Dry den and contained in MIL-STD-8785A 
Random gusts with the appropriate spectra are obtained by passing gaussian white noise 
through a linear filter. The filter' transfer functions for the three gust components 
are given in Table 9. The scale length (L^) and rms velocities (a.) have the following 
altitude dependence: 

1/3 

L = L = 145.0 h ' 
u v 


h > 533.4 m 
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Secondary Actuator Model 



Linear 


1 

1 

l < 
1 

First-Order 

Command j 

Second- 


Hysteresis 

Element 

J_ 

| 

Lag Plus 
Rate Limits 

■ 

1 

Order Lag 


! 

1 




j 

1 



Primary Actuator Model 


(Hysteresis 
Element 


,_J L_. 





Hysteresis element 


Figure 2. Simulation model of actuator. 
TABLE 9. GUST MODEL 


Gust Spectrum 

Transfer Function 


2L „ i 

G (s) 
ug 

u 1 

a u V L 

a , , u 

1+ y“ s 

a 


3L v 

l 1 + i r 18 

G is) 
vg' 

v a 

c v V L 2 

1+^ s 

Q. 


3L 

i . w 

L 1 + — 3 

G (s) 
wg 

w a 

w V L 2 

a - . w 

1 + 1T S 

a 
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L = L = Li = 1750 ft = ,533.4 m 

W U V 


1.25 


CT w 5, 25 " 1Og 10 3046 


h in meters 


This expression approximates the MIL-ST-8785A values from 30.48 to 18, 288 m 
For 0 4 h < 30.48, the value of h = 30.48 is used in the above equations. 


For design work, the transfer function for vertical gusts was approximated by a first- 
order filter: 



For the MLE concept, the vertical gust model incorporated in the Kalman filter was 
further approximated. Using the first-order differential equation model corresponding 
to (2), 



The low frequency content of the derivative was increased to yield 



SENSOR MODELS ! , • 

j 

Sensors were modeled as unity transfer functions plus additive noise. Sensor noise 
effects were modeled for the rate gyros, accelerometers, and servo position transducers. 
The noises were modeled as white noise shaped by a first-order low-pass filter. 

The filter time constants and rms values are given in Table 10. This spectrum 
approximates the bandwidth of most flight control sensors, as well as the nominal 
pre-filter used on sampled CAS signals. There is some difficulty in defining "sensor 
noise" because sensors in flight will measure unmodeled dynamics (body bending 
effects) and unmodeled disturbances (engine noise, etc.). The values presented here 
are considered representative of internal noise only. 


29 



TABLE 10. SENSOR NOISE PARAMETERS 


Sensor 

Si 

Correlation time 

rms value 

Pitch rate gyro 

0.01 sec 

0.0026 rad/sec 

Yaw* rate gyro 

' " - 0. 01 sec 

0.0026 rad/sec 

Roll rate gyro 

0.01 sec 

0.0131 rad/sec 

Normal accelerometer 

0.01 sec 

2 m/seo^ 

Lateral accelerometer 

0.01 sec 

0.03 m/!sec^ 

Servo position 

0.01 sec 

0.0007 rap 


Implemented as white noise shaped by a first-order lag with this 
time constant. 

DESIGN MODELS AND SIMULATIONS 

The limit cycle, model tracker, and maximum likelihood estimation designs were 
developed using several versions of the above models. Pitch and lateral-directional 
designs were done separately using uncoupled models. The pitch axis designs used 
only the short -period dynamics, elevator dynamics, and gusts. The lateral designs 
used a fourth-order aircraft model, aileron and rudder control effectors, and gusts. 
Other modeling details peculiar to a specific design will be covered in their specific 
sections of this report. 

A different simulation was used for the design of each adaptive concept. Each simula- 
tion used linear airplane equations of motion and constant parameters representing 
a specific flight condition. The simulations were tailored to facilitate design trade-offs 
for a specific concept. Additional details and rationale are provided in the sections 
devoted to the design details of each concept. 

The simulation used for design verification was the complete nonlinear six-degree-of- 
freedom simulation at NASA Langley Research Center. It operated in real-time on a 
CDC6600 and permitted man-in-the-loop simulation. Provision for driving actual 
F-8C servos with computer signals was included. Details of this simulation are 
contained in Reference 59. 
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SECTION 6 


FUNDAMENTAL LIMITATIONS OF ADAPTIVE 
CONTROLS FOR THE F- 8C 


Like most physical control processes, adaptive controllers for the F-8C aircraft cannot 
work arbitrarily well. There are fundamental limitations imposed by aircraft physics, 
measurement constraints, disturbances, and operational requirements which bound 
achievable performance. It is important to recognize and to evaluate these limitations 
before detailed adaptive designs are seriously undertaken. This establishes realistic 
"levels of expectation" for the performance of final designs and helps to define promising 
design configurations and initial design parameters for particular adaptive concepts. 

In the design program reported here, such fundamental limitations of adaptive controls 
for the F-8C were evaluated with separate identifiability studies. The purpose of these 
studies was to determine how accurately parameters in the aircrafts linearized equations 
of motion can be explicitly identified in an operational setting- -that is, with realistic 
sensor imperfections, realistic atmospheric disturbances, and, most importantly, 
with operationally tolerable test signal levels. Identification accuracy was selected 
as the key index of adaptive capability because it translates readily into control-loop 
performance for explicit adaptive concepts. It also bears a strong, though indirect, 
relationship to the performance of implicit concepts and can be computed theoretically 
without elaborate simulation evaluations. 

Separate identifiability studies were carried out for the pitch and lateral-directional 
axes. In each case, an identification model was first developed which included dominant 
short- period dynamics, disturbances, sensor characteristics, and trim conditions. 

This model was parameterized by expressing individual equation coefficients as functions 
of certain dominant aircraft parameters (q, M § , Mach) plus a set of small perturbation 
parameters. 

A linear system identification problem was then formulated to estimate these dominant 
and small perturbation parameters from on-line, closed-loop, input/output records. 

eo 

As shown in the literature, the theoretical accuracies attainable in such identifica- 
tion problems are bounded from below by the inverse "Fisher information matrix. " 
Hence, it remained only to compute this matrix for various operational situations and to 
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use its diagonal elements as fundamental identification accuracy limitations for the 
F-8C aircraft. 

Results show very limited identification potential. For example, out of six parameters 
describing pitch dynamics, only two can be readily identified. These are elevator 
effectiveness (M^) and a small perturbation 'parameter corresponding to pitching moment 
due to angle-of-attack (3VF ). The remaining parameters are known more accurately 
from a priori data than from on-line estimation. The estimation process does not 
significantly reduce their a priori uncertainties. Also, theoretical identification 
accuracies associated with and are no better than 10 to 20 percent depending on 
the length of data samples used. 

Similar results hold for the lateral- directional axes. In this case, only two out of 
16 parameters can be identified. These are dynamic pressure (q) and the perturbation 
parameter corresponding to yawing moment due to sideslip (N D ). On-line identifica- 
tion does not significantly improve a priori knowledge for any of the rest. Theoretical 
accuracies for the identifiable parameters are again in the 10 to 20 percent range. 

Tbe primary reason for these weak idehtifiability properties is the operational ground 

rule governing test signal magnitudes (Section 4). This ground rule makes explicit 

identification for adaptive control qualitatively different from, say, explicit identifica- 

42 43 

tion for post-flight data reduction. * In one case we are free to choose the magnitudes 
and shapes of control inputs so as to optimize identifiability. In the other case we must 
try to make the inputs as small as possible, hopefully to the point that they blend in and 
get lost among natural noise and gust induced signals of the control loop. This difference 
is fundamental to' the design work reported here. It affects the algorithms used for 
identification, their structure, the number of parameters identified/ and anticipated 
identification accuracies.- 

The remainder of this section provides ..support for the identification properties summarized 
above. Identification models and their parameterizations are presented for the pitch 
and lateral-directional axes, together with identification accuracy results as functions 
of data length* test signal magnitudes, and gust and sensor noise conditions. 
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PITCH AXIS IDEN TIFIA BILITY 
The Identification Model 

The following model was used for the pitch axis identifiability study: 




k .■ 0, 1, 2,.,. (t^j - t^) -A 

This is a simplified version of the aircraft models described in Section 5. It was 

derived by partitioning the 12-state models obtained from numerical linearizations into 

longitudinal and lateral- directional dynamics. The longitudinal states (q, id, u, 0. h, x ) 

were then transformed to (q, a * w/V, u, Y '«•». h, x) and partitioned into the fast 

short- period dynamics (q, a) and the slower states (u, Y» h, x). Coupling between these 

two groups was ignored. To complete the model, a first-order approximation of the 

Dryden vertical gust spectrum was appended; constant distuxbances (M q , g/V) were 

added to represent trim; and output equations were written for the gyro, q m » and 

the accelerometer, N . 

zm 
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the theoretical identification error levels reported in this section, provided that test 
signal levels are not changed. Conversely, they would permit the same error levels 
to-be achieved with proportionately smaller test signals. 


Model Parameterization 

The aoove identification model contains 11 parameters -- 

six coefficients describing pitch dynamics: 

M q , V , <Z a V>, M 6 , (Z 6 V> 

two gust level and trim terms: 


three initial conditions : 

q(t Q ). Wt 0 >+« g (t 0 >], « g (t o ) 

In the identifiability study, all 11 were treated as unknowns. However, because there 
is a good deal of a priori knowledge about their values and correlations, the parameters 
were not treated as independent, unrelated variables. Instead, they were modelled as 
dependent functions of two dominant aircraft parameters (M fi and Mach) plus a set of 
small perturbation parameters. The latter parameters (denoted by vector c) were then 
identified. 

The functional relationships between new variables and old were obtained by approximating 
scatter plots for the original parameters evaluated by 25 flight conditions given in 
Section 5. This gives very simple functions --linear (or quadratic at worst) in M fi with 
two sets of coefficients, one for subsonic and the other for supersonic flight. Scatter 
about these functi ons is small, providing tight a priori bounds for most of the new 

parameter values. The functions are summarized in Table 11, and their associated 
scatter plots are given in Appendix A, 


It should be noted that this parameterization is airplane-specific. The functions in 
Table 12 apply to the F-8C only. However, it is known that other aircraft also exhibit 

strong correlations between model coefficients, and, hence, similar tables can be 

readily constructed for other design applications. 
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TABLE 11. PITCH MODEL PARAMETERIZATION 


Original parameter 

New parameter 

New parameter a priori 
one sigma uncertainty 3 - 

M 

q 

= -0. 23 + 0. 028 Mg + Cj 
= -0.23 + 0.010 Mg + c. 

Mach < 1 
Mach > 1 

c i 

0. 065 

M 

a 

= (0. 61 + 0.92 c 2 ) Mg 

Mach < 1 

C 2 

0. 135 


= (1.53 +c 2 ) Mg 

Mach > 1 



V 

= (200 + Cg) ~Y~M§ 

Mach < 1 

C 3 

31.5 


= (260 + c 3 ) -V-Mg 

Mach > 1 



(Z«V) 

= (53 + c 4 ) Mg 

All Mach 

C 4 

5.0 

M 6 ’ 

= C 5 

All Mach 

C 5 

15.0 

(Z 8 V) 

= (7. 7 + Cg) Mg 

All Mach 

■ 

C 6 

0.750 

— - 


— — _ — 




a wg 

= 3.0 + c ? 

All Mach 

C 7 

1.50 

M 

o 

= 0. 017 Mg + c g 
= C 8 

Mach <1 
Mach > 1 

C 8 

0.125 

q^o). 

“ (t o ) + “g (t o > * “g (t o ) 


Same as original 

0.02 rad/sec, 0.02 rad, 
0 . 01 rad 


a Taken to be one-half of the + max uncertainty range in the figures in Appendix A. 




Accuracy Analysis Procedure 
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As reported in much recent literature, f ' one way to identify the vector £is to 
maximize the conditional probability density function of the observed outputs given the 
unknowns and the inputs, i. e. , 

6 = Arg p (q^, N zmfe ; k = 1, 2, ... . N | c = Ci & k » k = 0, 1, . . ., N-l)} 

^ * (15) 


Without repeating too much in the references, this maximization reduces to the 
following minimization problem: 

. 2 N 2 

c = Arg {™ in 1/2 [| | C-c 1 1 -p -1 + 2 <11 v. | | + l ndet B.)]} (16) 

C ~° F o k=0 B k 

where 

c = initial parameter estimate vector, 

— o 

P Q = initial error covariance matrix (Table 11), 

v k = residual (innovations) sequence of Kalman filter computed for parameters 
c ■ C, and 


B. = computed covariance matrix for the residual sequence. 


Let the term in brackets above be denoted by L(C, N), the log likelihood function. 
Then the required minimization operation can be accomplished iteratively using the 
following modified Newton- Raphson procedure: 



* These expressions assume high sample rates, so that B k “ R/At and VB^ = 0. 
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(17) 







Under suitable identifiability assumptions, this procedure converges to the true parameter 
vector, C 1 £.. The final iteration produces the second-partials matrix V^L (c., N) 
which is the so-called "Fisher Information Matrix. " Its inverse is a theoretical lower 
bound for the parameter error covariance matrix. Diagonal elements of the inverse 
are lower bounds on the variances of the individual parameter identification errors. * 

The identificability study basically consisted of carrying out the last iteration of the 
above algorithm for various noise, test signal, and data length conditions. The procedure 
is illustrated in Figure 3. For each case, the identification model (equations (5) and 
(6) with c: = c^) was used to generate a sequence of observations ^ = 1* 2, . . . , N}. 

This was done in closed-loop fashion, using a simple proportional-plus -integral pitch 
rate controller. The loop was excited by a small test signal <u.) and, for some cases, 
by a pseudorandom gust input (11). The sequence of observations was then used to 
generate the likelihood function, L(c t * N), its gradient, VL(c t , N), and its second- 
partials matrix, v 2 L(£ t< N). 



Figure 3. Accuracy analysis procedure. 


+ Strictly speaking, the actual identification error approaches the inverse from above 
asymptotically (as N gets large). Hence, the inverse represents a lower bound for 
the identification errors at all finite N. 

* In noise and/or gust-driven cases, L, vL, v 2 Li are stochastic quantities. For such- 
cases, these quantities were computed as sample averages of several (usually five) 
separate samples. 
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Three different types of test inputs were used: square waves, sine waves, and random 
signals. The first two were motivated by References 63 and 69 which indicate that 
cyclic test signals near the natural frequencies of the plant are good for identification 
accuracy. The last was chosen to provide a test input which is difficult to distinguish 
from atmospheric turbulence. Magnitudes for all signals were scaled to produce rms 
nominal acceleration levels in the range of 1 to 1.5 m/sec at the pilot station. These 
levels fall into the "comfortable” range of recent ride quality studies (for example. 
Reference 64) but are well above the so-called "threshold of detectability" which falls 
around 0.2 - 0,3 m/sec . 

The random test signal consisted of white noise shaped with the following linear filter: 

G(s) '= ~ g 2 < 18 > 

s + 2£uj s + u) 

The frequency (u>) was selected near the short period of the aircraft; the damping ratio 

2 

C was set to 1. 5; and k was adjusted to provide the desired 1.0 - 1.5 m/sec acceleration 
at the pilot station. 

The study also addressed reduced-parameter identification, i. e. , identifying fewer than 
11 error parameters. This is an important practical question since all 11 parameters 
comprise a fairly high order identification problem whose solution is probably not 

feasible on board. For reduced- parameter identification, the error covariance matrix, 

2-1 

(7 L) , must be modified slightly to account for the errors produced by the presence 
of incorrect unidentified parameters. The modified matrix is derived as described in 
the following paragraph . 

T 

Let the parameter vector £be partitioned into two groups, [c^, c^l • The first group 

is to be identified, while the second is to be neglected. Then the likelihood function 

T 

near the true parameters, c = [£ tl * , can be written as 

LCC, N) = Me N) + i II .{ - c II 2 +6 (19)* 

-t 1 -t v 2 L 

VL = V 2 L (C - c ,) + e 


The first-order term vanishes because L(g, N) is minimized at Q = c^. 
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where 6 and e represent higher order term and random errors. 

A reduced parameter Newton- Raphs on minimization procedure will now iterate on 
g* according to Algorithm (17) until the gradient components vLj (not the whole vector 
VL) is zero. This condition is satisfied when 

- £ n = - (V 2 L n ) _1 [V 2 L 12 (G 2 - c t2 ) + 6l ] <2D 

Hence, the error covariance of the converged estimate will be 

Co. <£j> i EW 1 , - o u ) (c\ - £„) T ] ■ . 

. (. 2 L)"J [(7 2 L) 11 + (. 2 L) 12 P 22 <7 2 L) 21 ] (V 2 ^ 

where P 22 = B[(C 2 - £ t2 > <C 2 - £ t2 ) T ] is the a priori covariance matrix of the unidentified 

parameter, components and (V 2 L) U can be shown to be the covariance matrix of * r 
Equation (22) agrees with the full parameter identification theory if P 22 = 0. 


A naly sis Results 


Results of the pitch axis accuracy study are summarized in Table 12 and Figures 4 through 
7. The principal result is Table 12 which shows before and after accuracy comparisons 
for a 0. 13 g square wave test signal and true parameters corresponding to Flight 
Condition 1. The table compares a priori parameter uncertainties from Table 11 
(the "Before" column) with parameter uncertainties obtained after identifying with five 
seconds and ten seconds of observed on-line data (the "After" columns). Scanning these 
i numbers shows that only c 2 . c & , c g , Cg, and c 1Q can be improved by more than a factor 
of two. Among the first six dynamic parameters, only c 5 (M fi ) is improved dramatically, 
c (small perturbation for M # ) improves marginally, and the rest are barely touched 
at all. The corresponding percentage accuracy at ten seconds for M 6 is 
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TABLE 12. BASIC IDENTIFICATION ACCURACY 

1/ 

0. 13 g square wave test signal 
for Flight Condition 1 
40 samples /second 


Parameter 

c i 

(M q ) 

C 2 

(M «> 

C 3 

(V) 

C 4 

(z a v) 

C 5 

(M 6 ) 

C 6 

(Z 6 V) 

c 7 

(cr ) 
wg' 

C 8 

< M o> 

C 9 

<v 

C 10 


C 11 

<« > 
go 


After 


One sigma uncertainty (t = data length) 


Before 


sec 


0.065 0.0630 0.0620 

0.135 0.0650 0.0560 


0.750 


0.125 


= 5 sec 


0.750 


0.0063 


0.0066 


0.0058 


0.0063 


t = 10 sec 


0.0620 

0.0560 


1.46 

0.740 


0.0050 


0.0067 


0.0053 


0.0061 













M 


fit' 


1.46 

14.8 


10 percent 


(M Ci = true value of M.) 
o.t o 


and for M it is 
a ■ 


'5 


M 


6t' 


0.61 


S c 2 


6t 


) (q - gj) Cl 10 percent 


If we look at the parameter alone, its identification accuracy depends on data length 
(as seen above) and also on the true parameter vector c^. This is shown in Figure 4 
which presents results for three flight conditions plotted as a function of data length. 

The most significant point to notice is that M g errors are more or less proportional to 

the true M values. This is consistent with the figures in Appendix A which show that 

8 

a priori uncertainties are also proportional to the true values. , It permits us to express 
subsequent identification errors as normalized quantities. The data length dependence 

of M errors is also worth noting. Errors are reduced to 20 to 30 percent almost 

8 

immediately and then improve very slowly with time from there on. As shown later in 
Section 8, it is not possible to use more than five to ten seconds for identification because 
the aircraft parameters can change significantly over such time spans. Hence, the 
accuracies in Table 12 and Figure 4 represent practical lower bounds on explicit on-line 


identification performance. 


The possibility of reduced-parameter identification is examined in Figure 5. This figure 
was obtained by selectively deleting error parameters from the list of identified unknowns 
but accounting for their uncertainty in subsequent error calculations. The results 
show that many parameters are insignificant with respect to accurate identification. 

In fact, only itself, , and the initial condition (a + a g^ 0 rea ^y matter much to 
the accuracy of M g . This is an important conclusion because it justifies simplified 
identification procedures for use onboard the F-8C. 

* 

The results presented so far include sensor noise and data uncertainties, but not gusts. 
The latter are added in Figure 6 and improve identification slightly. Figure 6 also shows 


* It should be noted that the maximum likelihood procedure was in all cases set up to 
expect gusts even when none were included in the true aircraft responses. 
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0.13 g's square wave test signal 
for FC1 , 10 and 17 



Data length 


Figure 4. Identification accuracy for M 


Parameter identified 
c 5 c 8 c 10 c 7 c 2 C 4 C 1 C 3 

M (a+a a „ M Z V M V 
' 6 "0 ' g 0 oig a a q 

XX XXX XXX 

XXX XXX 

X X X X X 

XX X X 

XX X 

XX 

X 

g* Figure S. 




Normalized identification error 







effects of test signal magnitude. As expected, larger test signals improve accuracy, 

at least up to a point. About 10 percent error remains even with very large test inputs. 

This is caused by errors in the unidentified model parameters. (Note that Figure 6 

uses reduced-parameter identification.) For purely cyclic test inputs, these parameters 

have identical effects on M. error over each period of the input, and, hence, no improve- 

o 

ment occurs with time. 

The latter effect can be reduced by gust inputs and by using test signals which are not 
cyclic. This is shown in Figure 7, where random test signals are "compared with square 
waves and sine waves. Both cyclic signals are seen to be inferior to random commands. 
Compared against each other, they are roughly equivalent. 



Figure 7. Effects of test signal type. 


LATERAL- DIRECTIONAL IDENTIFIABILITY 

Similar accuracy analyses were also conducted for the F-8C lateral-directional axes. 
Results for this case are reported here in abbreviated form. They point to the same 
basic observations— only a few parameters are readily identifiable on-line and those 
have only limited accuracy. 
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The identification model for the lateral- directional axes is shown in Figure 8. Like 
pitch, this model was developed from the original 12-state aircraft linearizations in 
Section 5.; These were first partitioned to separate lateral-directional dynamics* 

The lateral states (p, r, v, 0, i|r, y) were then transformed to (p, r, 3 = ~ , y), 

and the fast modes (p, r, g) were isolated. 

Finally, the model was appended with an approximate lateral gust model, a disturbance 
term (g) to represent geos 9 sin0 gravity inputs, and algebraic sampled output equations 
for the roll and yaw rate gyros and for the lateral accelerometer. 

Like the vertical, spectrum, the lateral gust spectrum was approximated from the complete 
Dryden model in Section 5. The approximation follows Equations (7) through (10) exactly 
(except, of course, for symbols) but adds the following additional simplification: 

ar V*’ ' -<ir> V> + > w‘> 

(23) 

- ,,<i) 

In effect, this added simplification turns the gust derivative into white noise (i.e. , a 
flat spectrum). This is pessimistic in the sense that the actual derivative spectrum 
is flat at high frequencies only and falls off below (V/L^) rad/sec. 

The gravity disturbance (g) was added in lieu of the bank angle state which was not 
considered to be available for the study (Section 4). It was modelled as a first-order 
random process with bandwidth a = 1 sec and rms level = 9.8 m/sec^. The sensors 

g 

were again assumed to have wide-band noise outputs so that low frequency sampling 
will produce white discrete noise sequences. The rms sensor noise levels were identical 
to the pitch study. 

The lateral-directional model was parameterized in terms of two dominant aircraft 
parameters--q and Mach--and a new set of small perturbation parameters. The trans- 
formations are again simple linear functions (or quadratic, at worst) with different 
coefficients for subsonic and supersonic flight. The functions are summarized in Table 13. 


48 






CO 


Figure 8. Lateral- directional identification model. 



I u* 


TABLE 13. LATERAL-DIRECTIONAL MODEL PARAMETERIZATION 


Original parameter 

New parameter 

New parameter 
one sigma a 
uncertainty 

q 


c i 

Ail Mach 

■ 

c i 

0.20 c 2 . 

V 

= 

(41 + c 2 ) ~\fq 

Mach < 1 

C 2 

7.0 



(54 + c 2 ) Vq" 

Mach > 1 



L 

P 

= 

-1*, 5 - 0. 0071 q + c 3 

Mach < 1 

C 3 

0.60 



0.5 - 0.0071 q+c 3 

Mach > 1 


0.30 

L 

r 

= 

0. 44 + c^ 

Mach < 1 

C 4 

0.26 



-0*03 + c^ 

Mach > 1 


0.020 

L |3 

= 

-10 - 0.12 q+c 5 

All Mach 

C 5 

5.0 

N 

P 

= 

-0.090 - 0.00012 q+c- 

D 

Mach < 1 

C 6 

0.025 



-0.025 - 0.00012 q+c 6 

Mach > 1 


0.009 

N r 


-0*10 - 0.00090 q+c^ 

Mach < 1 

C 7 

0. 040 



0.02 - 0.00090 q+c ? 

Mach > 1 


0. 030 

N |* 

- 

-2* 3+0. 019 q+c Q 

o 

All Mach 

C 8 

1.20 

V 


(-0. 52 + c g ) q 

All Mach 

C 9 

0.060 

L 6a 

= 

(0.075 + c 1£) ) q 

Mach < 1 

C 10 

0.008 



(0.032 + Cjq) q 

Mach > 1 


0. 004 

, 

L 5r 

= 

(0.029 + j) ^ 

Mach < 1 

C 11 

0.004 



(0. 014 + c 11 ) q 

Mach > 1 


0.002 

N 6a 


(0.0038 + c 12 ) q 

Mach < 1 

°12 

0. 0006 



(0.0013 + c 12 ) q 

Mach > 1 


0.0004 


' 'cl 

Taken to be one-half of the + max uncertainty range of the scatter plots* 
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TABLE 13. Concluded. 


Original parameter 

New parameter 

New parameter 
one sigma a 
uncertainty 

N 

or 

(-0.013 + c 13 ) q 

Mach < 1 

C 13 

0.0020 


(-0.005 + c 13 ) q 

Mach > 1 


0.0004 

Y 6 a V = 

(0.0094 + c 14 ) q 

All Mach 

C 14 

0.0028 

Y 8r V = 

(0.110 + c 15 ) q 

Mach < 1 

C 15 

0.016 


(0.042 + c 15 > q 

Mach > 1 


0.004 

P< t o ) = 

1 

1 50 
1. cT 

All Mach 

C 16 

0.017 

P(t 0 ) * 

C 17 

All Mach 

C 17 

0.017 

Mo> + W C 18 

All Mach 

C 18 

0.017 

2 ( V = C 19 


All Mach 

°19 

32.2 


a Taken to toe one-half of the + max uncertainty range of the scatter plots. 
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The primary accuracy results for the lateral-directional axes are summarized in Table 14. 
This table provides a before and after comparison of a priori parameter uncertainties 
(t - 0) with uncertainties obtained after identification with two, five, and 10 seconds of 
on-line data. As in the pitch case, only a few parameters are readily identifiable- -two 
out of 16 dynamic parameters (q and c g (N p )) and three of four initial conditions. 

Accuracies are again in the 10 percent range. These results were obtained for closed- 
loop operation with a simple roll rate control loop. The loop was excited by random 
test signals generating 0. 05 g's rms lateral acceleration. The lower level of lateral 
acceleration is: consistent with ride quality requirements. 64 

F-8C IDENTIFIABILITY SUMMARY 

The pitch and lateral-directional identifiability studies led to four major observations; 

• Only a few F-8C parameters can be realistically identified in an operational 
on-board setting. These include M g and M ft in pitch and q and N g in lateral. 

• Accuracies no better than 10 to 20 percent can be anticipated for the identifiable 
parameters. This accuracy range applies only for the sensor noise, gust 
levels, and test signal magnitudes assumed in the study, of course. It can be 
approximately scaled to other assumptions, however, by noting, in Algorithm 
(9), that the covariance matrix is directly proportional to B k (sensor noise 
variance) and inversely proportional to Vv w T (signal level squared). Using 
this scaling, for example, we can expect to get the same accuracy with lower 
test inputs if we have proportionately cleaner sensors. 

Reduced-parameter identification can successfully find identifiable parameters 

with little loss of accuracy. The parameter set in pitch can be reduced to as 

few as three parameters --M g and two initial conditions, M and (a + a )— 

and still provide good M. estimates. * 

0 

Random test inputs appear to be superior to cyclic inputs for reduced-parameter 
identification. 

These results have significant impact on the adaptive control concepts formulated in 
later sections. 
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TABLE 14. LATERAL AXIS BASIC IDENTIFICATION ACCURACY 



One sigma uncertainty (t = data length) 


5 sec 



60. 

7.0, 

0.60 

0.26 

5.0 

0.025 

0.040 

1.20 

0.060 

0.008 

0.004 

0.0006 

0.001 

0.0028 

0.016 


0.017 

0.017 

0.017 

32. 


33. 

6.7 

0.58 

0.26 

4.1 

0.025 

0.040' 

0.66 

0.060 

0.008 

0.004 

0.0006 

0.0009 

0.0028 

0.016 


0.011 

0.003 

0.004 

18. 


26. 

6.6 

0.57 

0.26 

4.1 

0.025 

0.038 

0.61 

0.060 

0.008 

0.004 

0.0006 

0.0009 

0.0028 

0.016 


0.011 

0.003 

0.004 

18. 


23. 

6.3 
0.55 
0,26 
3.7 
0.025. 
0.038 
0.58 
0.060 
0.008 
0.004 
0.0006 
0. 0009 
0. 0028 
0.016 


0.011 

0.003 

0.004 

18. 


0.052g rms random test signal 
0. 012 rad/sec rms roll rate 
0. 005 rad/sec rms yaw rate 










SECTION 7 


THE CONTROL LAW STRUCTURE 


The optimal control laws of the CCV package were used as a basis in designing the 
adaptive control laws. The major difference being that the gains scheduled with air- 
data quantities are replaced with gains determined by the adaptive concepts. 

The pitch and lateral- directional control laws will be described below. These designs 
are explained in greater detail in Reference 2. 

PITCH CONTROL AUGMENTATION SYSTEM 

The structure of the pitch control law is shown in Figure 9. It is a C*™model~following 
design that was designed with linear quadratic optimal control theory. Signal shaping 
and an integrator were included in the optimal control problem formulation to satisfy 
design specifications. The optimal design has been redrawn here to illustrate the 
sensors used and the signal shaping employed. The C*~ quantity is formed by a fixed 



Figure 9. Pitch CAS structure. 
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sum of normal acceleration and pitch rate multiplied by (the "crossover" 
velocity), a constant of 324 ft/sec. The normal acceleration is lagged and a lead- 
lag compensation is applied to pitch rate. Values for these filters resulted from 
the quadratic design. Apparent netural speed stability is provided by generating 
an integration in the forward path of the elevator control system. This is 
mechanized by positive feedback of the elevator secondary servo position. In 
actual operation, the average position is used since the F8-C has a split horizontal 
tail surface. The loop gain (G^ + ) was previously scheduled with dynamic pressure 
as shown in Figure 10. For the two explicit adaptive concepts, this gain is adjusted 
from estimates of elevator surface effectiveness. For the implicit concept 
(high-gain limit cycle), the gain is adjusted to remain at or near critical gain. 

The specific values of the filters are given in Figure 11 for the digital pitch CAS 
at a 32 sps sample rate, which was the rate on the F-8C simulator. The flight 
computer will operate at 50 sps. The analog equivalents of the digital filters 
are indicated in parenthesis. 

LATERAL-DIRECTIONAL CONTROL AUGMENTATION SYSTEM 

The structure of the lateral-directional CAS used with each of the adaptive 
concepts is shown in Figure 12 in a standard block diagram form. This design 
resulted from a reduced measurement solution of an optimal control problem. 

The full-state optimal controller was designed to follow an explicit roll rate 
model and to minimize sideslip and lateral acceleration for roll stick commands. 
The resulting design produced good turn coordination and Dutch roll damping over 
the entire flight envelope. 

For implementation on the F-8C, the full-state design was simplified to use a 
reduced set of measurements. The reduced measurements consisted of p, r, $ a » 
and N . For the CCV package, the feedback and feedforward gain schedules 
were modified to be a . simple function of angle-of-attack. This eliminated any 
need for true airspeed or attitude measurements. Details are contained in 
Reference 2. -~~ 

For use with the adaptive control laws, the feedback gains of the reduced measure- 
ment CAS were grouped into a single-loop gain that is now a. function of surface 
effectiveness (or dynamic pressure). (Refer to Figure 13. ) 
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0.1743 + 0.1743Z 
1 - 0.651 33Z' 1 


0.01 56+0. 01 56Z~ 



0. 03587+0. 03587Z' 1 
1 - 0.9283Z” 1 . 


0.3639 


1 


(q.42S + 


0.0156 + 0.0156Z 

TTTT— 


Elevator servo 
-*► command (rad) 


0.03587 + 0.03587Z 
1 - 0.9283Z"* 
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The two explicit adaptive concepts determine the loop gain from surface effective' 
ness estimates in the pitch axis. The limit cycle operates with a gain changer in 
the roll axis that keeps G LATi near critical gain. The yaw axis loop gain is 
slaved to the roll axis. The digital lateral-directional design is shown in 
Figure 14 for a 32 sps sample rate. 
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Figure 14* Lateral-directional CAS, 
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SECTION 8 
THE MLE DESIGN 

One of the three candidate adaptive concepts selected at the outset of the design program 
was intended to explore modern estimation-theoretic methods for explicit parameter 
identification. These were to be combined with simple algebraic gain update calculations 
to make up the complete adaptive concept. After some initial experimentation with 
two identification approaches, maximum likelihood estimation and nonlinear filtering 
(the latter largely unsuccessful), it was decided to base the concept on maximum 
likelihood methods. The resulting adaptive controller design is described in this section. 

MAXIMUM LIKELIHOOD IDENTIFICATION 

We have already introduced the general process of maximum likelihood identification 
in Section 6. The principle is to find parameter estimates which maximize the a 
posteriori probability distribution for the observed outputs conditioned on the unknowns 
and the measured inputs, i. e. , 

£ = Arg [max pfyj, y 2 » yjgl — = P U q* u 1* *** U N-1^ 

When the unknowns are constant and the plant dynamics are linear, this maximization 
problem leads to the solution shown in figure 15. # The solution consists of a 

Kalman filter designed for the true system structure but with parameters equal to- an 
estimate c = C. The filter tracks the true system outputs and generates a residual 
| sequence = y fc - y fc k = 1, 2, ...}. This sequence is accumulated into a 

likelihood function, 

L(c, N) i - to p (y 1 , y 2 , ... , y N j £ = Cs V u l* * “n-1* (25) 

! which is then minimized over the parameter estimate £. 

At first glance, this solution appears ideal, for on-board applications. The Kalman 
I filter is relatively simple and runs recursively, processing data samples as they appear. 
I The same is true for the likelihood accumulation operation. The difficulty, of course, 
is the last step of the solution: the minimization. This requires repeated or parallel 
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. Inputs 



Figure 15. Maximum likelihood estimation. 

processing through the data and adds significantly to computational complexity. 

Approaches to minimization and the equations involved are discussed further below. 

Minimization Algorithms 

Two algorithms were considered for the likelihood minimization operation: 1) iterative 
Newton- Raphson calculations, and 2) parallel non-iterative calculations. The first 
algorithm is illustrated in Figure 16. It begins by collecting a sequence of input/output 
data, 

Y N = {y l' y 2* * y N^ 

(26) 

U N = U l* **• * °N-1^ 

The Kalman filter is then run with parameters c equal to a priori estimates, £ = C° £ 0 * 
This generates the likelihood function, L. At the same time, a set of sensitivity 
equations are processed which permit calculation of first- and second-partial derivatives. 
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Figure 16. Newton-Raphson minimization (Algorithm 1). 


2 

VL and V L (some form of approximation is usually used to simplify equations for the 

latter). These derivatives are then used to obtain a new parameter estimate using a 

standard Newton-Raphson formula. The filtering, accumulation, and Newton-Raphson 

operations are performed repeatedly for the same data set until convergence is achieved. 

The data set itself is usually kept current by a "sliding window" process which increments 

the lower and upper indices of the input/output sequences (Figure 16) while keeping 

the total number of samples, N, unchanged. The overall algorithm has been implemented 

42, 43 

successfully for various post-flight data processing applications . 

The second algorithm is illustrated in Figure 17. It achieves the same end result 
as Algorithm 1 but replaces iterative calculations with parallel ones. The sequence 
of input/output observations is sent simultaneously to M Kalman filter channels, each 
with its own sensitivity calculations and likelihood accumulations. The channels are 
distinguished by their assumed parameter values. Each one operates with a different 
parameter estimate, c - and hence computes, L^, and at a different 
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Figure 17. Non-iterative parallel minimization (Algorithm 2) 













fixed point in parameter space. The likelihood functions at these points are then 
compared to find the approximate minimum point, and a single Newton-Raphson step 
is taken from there to estimate, the true extremum. 

As in Algorithm 1, some procedure must generally be added to keep the processed 
data current. This is done in Figure 17 by high-pass operations which provide 
exponential de-weighting of past accumulated data samples. This can be contrasted 
with the above-mentioned sliding window process, which provides uniform weighting 
for the past N samples and zero weighting for all earlier ones. While not proposed 
in the exact format of Figure 17, general parallel estimation structures have been 
suggested in several references. ^ 

Assuming that we have a well-posed identification problem (i. e. , L(£, N) has a unique 
minimum at £ = c^ as N -> ®), it is clear that both algorithms will perform the desired 
estimation function. They have important differences, however, which favor Algorithm 2 
for on-board application. The chief differences are: 

• Recursi venes s - - Algorithm 2 is completely recursive. It processes data 
samples as they appear and needs no data base storage and data base manage- 
ment operations as required in Algorithm 1. 

• Iteration- -Algorithm 1 is iterative. This leaves both the data processing 
flow and computing time uncertain, and it raises convergence issues. In 
contrast, the computation flow and time of Algorithm 2 are completely defined 
and inviolate from cycle to cycle. 

0 Accuracy- -Both algorithms can find the minimum with arbitrary accuracy 
if given enough iterations or enough channels, respectively. However, for 
equal computing time (the number of iterations equal to the number of channels). 
Algorithm 1 is probably better. 

As criteria for on-board application. Items 1 and 2 far outweigh the importance of 
Item 3. This is especially true because the location of the minimum of a realization 
of L(£, N) is itself so inaccurate (Section 6) that errors made in finding the minimum 
contribute little to overall identification error. Hence, it was decided to proceed with 

the parallel MLE algorithm, 

* The enclosed superscript (i) is used to denote variables and parameter values for 
Channel i. 
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Filter , Sensitivity, and Likelihood Equations 


The basic data processing equations which must be solved by each filter channel are 
summarized in this subsection. They are stated in terms of general symbols correspond- 
ing to the following discretized plant equations: 

x k+i = Ax k + Bu k + r11 k 

(27) 

y k = Hx k + Du k + N ?k 


The matrices A, B, F, H, D, and N should all be thought of as dependent functions of 
the parameter vector £. Then the channel equations are: 

67 

Filter equations, 

x k+i =Ax k + B \ 


\ + i = y k+i - ^k+i - ^k+i 
x k+i = x k+i + K k+i \+i 


(28) 


Sensitivity equations for each component c of c (derived by differentiating (28) 

P 

with respect to c^), 

V x. = Av x. + (V A)x. + (v B)u. 
p k+1 p k p k p k 


V v, , = -HV x. - - (v H)x, Ll - (v D)u. ... 

p k+1 p k+1 p k+1 p k+1 


(29) 


V x. , . = V x. , . + K. , i V v. , . + (V K. , . )v. . . 
p k+1 p k+1 k+1 p k+1 p k+1 k+1 


Filter gains. 
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p, J . 1 = a p. a + r r 

k+1 k 


B k« ■ ® p k« « T + n nT > 


— T -1 
K, _ L1 = P, H B. 
k+1 k+1 k+1 


( 30 ) 


P k+1 = P k+1 " K k+1 H P k+1 
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Likelihood accumulation. 


] W^ + 1/2 [ i | a + iI 


+ in det B 


B 


-1 

k+1 


k+l J 


(31) 


VL k+ i = „vL k + W k+1 T B k+1 _1 v k+1 + 1/2 Trace (v^ \ +1 - B k+1 ) V(B k+1 > 


y2 L k+1 * = ^k + Vv k+1 T B k+l“ 1 V \+l 


with jjl - exp (-AI/t) for exponential de-weighting of past data* The choice of t is 
discussed later in this section* 


These equations warrant two explanatory comments. First, there are no sensitivity 
equations for the filter gains. This is because the matrix VK was computed by numerical 
differentiation throughout the design program, i. e, , 

VpK = [K(c + X p e p ) - K(c)]/\ p 

where e is a unit vector in the i-th coordinate direction and was chosen small 
P ■ 

compared to the range of c p . 

Second, since V 2 L must be a derivative of VL, it is evident that the last equation in 
(31) includes approximations. All second-partial derivatives and products of derivatives 
have been ignored. This is a common approximation for so-called modified Newton- 
Raphson procedures 2 ^ and has the important advantage of eliminating second-partial 
derivative sensitivity equations. It is also common practice to eliminate the last term 
of the VL equation and to replace the filter gain equations with their steady state solutions. 
These additional approximations are considered later. 

IDENTIFIER DESIGN FOR THE F- 8C 

In view of the identifiability results in Section 6 and the known scheduling requirements 
of the control laws (Section 7), a specific implementation of Algorithm 2 was developed 
for the pitch axis only and for a reduced-parameter set with three components. Design 
details and trade-offs for this identifier are presented below. Key design issues include: 
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Identification models 


• Channel selection, 

• Kalman filter design, 

• Adaptation to noise statistics, and 

• Likelihood filters. 

Identification Models 

Following Section 6, a simplified short-period model was chosen to represent the pitch 
plant dynamics (27). 


State equations: 
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This model includes short-period dynamics (q, or), turbulence trim states 

(g/V, M q ), and noisy sampled algebraic sensors for pitch rate, normal acceleration, 
and elevator surface position. Turbulence is again modelled as an approximated 
Dryden spectrum without the low frequency drop-off in its derivative, i. e, , 


_d_ 

dt 





w 


a + a 
g w 




Tl 


Hence 



(34) 



The trim terms (which according to Section 6 must be included for accurate reduced- 


parameter identification) are appended as slowly varying states with Brownian motion 
representations. The term g/V was given a variance growth rate proportional to roll 
angle variance, and M q has a variance growth rate proportional to the variance of the 
trim pitching moment due to control. Hence, 



a 


9 



(35) 



trm 


M 6 



All sensor noises are modelled as independent and identically distributed random 

variables from sample to sample. Their magnitudes cr q , a nz , <7 g as well as the 

disturbance magnitudes o , cr , a were treated as design parameters. 

w 9 °trm 

Using Section 6, Table 11, the identification model was parameterized in terms of M § 
and two small-perturbation parameters for and V. The parameter M g was chosen 
because it is most identifiable and also most important for control. was chosen 
because it falls next in the hierarchy of identifiability. Velocity was included because it 
is desirable for lateral control. However, since V is not readily identifiable from the 
pitch fl-Hs alone, we cannot expect much in the way of accuracy. The rest of the 
perturbations in Table 11 were assumed to be zero for reduced-parameter identification. 
The parameterization was also modified to account for quasi-static flexibility of the 
airframe (not included in Section 6). This phenomenon reduces all control effectiveness 
terms at high dynamic pressure (high M g ) and can be accounted for by expressing the 
model in terms of an "un-flexed" M go parameter instead Of M g directly. The resulting 
parameterization is summarized on the following page. 
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M 6 


M. (1 +0.016 M +0.0002 M, 2 ) 
5o 80 6o 

i r 

M q 


! -0.23 + (0.028 - 0.017 c 0 ) M c 

! Z 0O 

M « 


; (0.61 +0.92 c„) M s 

£ OO 

V 


(200 + c 3 ) M 6o 

Z « V 


53 M 

6o 

Z 5 V 


7.7 Mg 

W _ 




c 



f 0 subsonic 
c *** \ 

2 1 supersonic 


C 3~ 


0 subsonic 
60 supersonic 


( 36 ) 


(37) 


Channel Selection 

The three-dimensional parameterization above generates a parameter space illustrated 
in Figure 18'. The space Consists of two regions, one for subsonic and the other for 
supersonic flight. There is also a loose interconnection through which the aircraft 
migrates rapidly in the transonic flight regime. The problem of channel selection is to 
choose both the number and the location of points in this space at which to operate 
Kalman filter ajid sensitivity equations. Criteria for this choice are derived from the 
interpolation function of the channels. We obviously want as few channels as possible, 
but they must be close enough together to be able to interpolate throughout the space. 


For the F-8C problem, a bit of experimentation showed that the c and c. directions 

a *5 

could be readily interpolated from channels located at their nominal values, i. e. , 

(0, 0) or (1, 60). The direction, however, requires several channels, each 

positioned so that it interpolates over an region which is a fixed percentage of 

its nominal value. This follows from results in Section 6 which indicate that accuracy, 

and hence curvature of the likelihood function, is proportional to true M e values 

o 
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Figure 18. MLE channel selection, 



(percentage accuracy is constant). Using these findings and a 50 percent interpolation 

range for each channel gives the following candidate M values : 

oo 


mJ 4) . -40/1.5 

= -26.7 

“so 6 ’ ’ M So <4>/2 - 25 

= -11.9 

M 5 o <2) ‘ “ jo ' 3 '/ 2 - 25 

= -5.27 

M 6q (1> ■ Mj 0 (2> /2,25 

= -2. 34 


These values must be combined with (c^ c 2 > = (0. 0) for subsonic flight and (c^ c 2 ) = 
(1* 60) for supersonic. A total of five channels are required — one for supersonic, 
which has an M gQ range entirely covered by interpolation about M go ^ = M^ 4 *, and 
four for subsonic. The resulting channel locations are shown in Figure 18, 

Kalman Filter Design 

According to Algorithm 2, a fixed set of Kalman filter, and sensitivity equations must 
be operated at each set of parameter values shown in Figure 18. While the design of 
these filters is straightforward, the following features of the F-8 design are worth 
noting. 

Time Invariance — A few experiments showed that steady state filter gains achieve 
almost the same speed of convergence (of the identifier) as time-varying gains. Hence, 
Equations (30) were replaced by their steady state equivalents (set k+1 = k) and solved 
off-line. Stored values of K and yK are used in the on-line identifier. 

Self ^Initialization- - When the adaptive mode is engaged, the filter and sensitivity 
equations in each channel are initialized on the basis of current sensor outputs 
(<l m » n zm » ® m )« This step avoids long channel transients associated with the slow 
trim states, g/V and M q , The initialization assumes that the aircraft is in equilibrium 
at the time of engagement (i.e. , x = 0), which gives the following equations for each 
channel: 
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q(0) 


V 


®(0) 

= 






(38) 

v” (o> 

= 

-' A 21 Sn + (A 22 ■ *> “ <0) +B 21 S m 1/A 23 

'• " • ( 

M (0) 

O 

S 

V * A 12 * <0 > + B 1X ‘ J ,A U 


Vq(0) 

= 

0 


va(0) 


-[(Z a V) «(0) + V(Z 6 V) 6 m ]/(Z a V) 




(39) 

v|(0) 

= 

-[VA 21 q(0) + 7A 22 Of(0) + (A 22 -l) v»(0) + VB 21 5 m ]/A 23 


VM (0) 
o 

= 

-[vAjj q(0) + VA 12 a(0) + A J2 7a(0) + VBjj 6 m ]/ A j 4 



Noise Statistics - -The filter gains for each channel were designed for the following 
assumed noise statistics: 

Sensor Noise Disturbances 

a = 0*0026 r/s = 6. 6 ft/s (1* 83 m/s) 

a = 0.02 g's (0.2 m/s 2 ) cr = 0.017 rad 

nz 

as trm ~ rad 

f ■ .. . . ■■■■■'■■ . . 

<r = 0,0008 rad 
6 


The Ain symbols in these equations refer to elements of the discretized dynamics 


matrix corresponding to Equation (32)* 

The Kalman filter Equations (28) treat u = 6. as perfectly known. The actual measure- 
ment, of course, is noisy. This was recognized by adding the control measurement 
error as an extra random disturbance in the discrete model Equation (27), i.e.. 


k+1 


= Ax k + Bu mk + (r \- B? 6k ) 


(27'). 
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Since these numbers are fixed for all operating conditions (i. e. , with and without sensor 
noise, with and without gusts, with slow or fast trim changes, etc, ), they represent 
experimental compromises over many simulation runs. As seen later in the performance 
summaries, some performance benefit can be gained by adapting these statistics to the 
environmental situation. 


State Estimate Interpolation - - Since each filter operates at a fixed location in parameter 
space, no one of them provides proper state estimates unless the parameters happen 
to fall close to one of the channels. This difficulty is alleviated, at least for the & 
estimate, by interpolating between adjacent channels. Because trim angles -of -attack 
are inversely related to M gQ , the following interpolating function was chosen: 


S (M 6o > - a, 



(40) 


The coefficients and a 2 for this function are computed from estimates of a at two 
adjacent channels (i, j) and then inserted back into (40). The resulting formula is: 


• (M 5o- *> 


<" t0 • “ <“*• «> + (M «o 1 


A A 

m ■ > "(Mr;, t) 

6o 5o 


< M S „ U> - 


M U) ) 
6o 


(41) 


Adaptation to Proportional Noise Statistics 

We noted in the filter design section that fixed statistics were used to compute filter 

gains. This is desirable because it generates an invariant set of gains for each channel 

which can be computed off-line and stored for on-line use. Invariant gains are actually 

obtained under slightly less restrictive circumstances — namely, when the disturbance 

and noise statistics remain in constant relationship to one another. This means that 

a filter for disturbances statistics (a , a , a , a r ) and sensor noise statistics 

w cp 6trm 6 

(<jq, a nz ) remains unchanged when those statistics are changed to (aa^ oc^ p aa 6t rm * 
acr g ) and (cro^, crCr nz ) for arbitrary a. As a result, the total identification algorithm 
designed for one set of statistics can still be valid when those statistics are scaled 
up or down, provided only that the likelihood functions are also scaled. This can be 
done adaptively as described in the following paragraphs. 

Assume that alii statistics are known to within a constant scale factor a. Then the 
likelihood function, in terms of the usual unknowns (£) and the new unknown a, is given by 
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N -1 

L = 1/2 2 (v^ T ~~2~ \ + An det cr 2 B) 

k=l a 


( 42 ) 


We find estimates for c and a in the usual way by letting 


0 = vL = 2 [? v k T ~~ 2 ~ \ + !/ 2 Trace (v fe v k T -a 2 B) Vf^-jj-)] 

k=l a 


,-l 

T 


N "1 

0 = |^ = 2 1/2 Trace (v, v. T -a 2 B) (-2 ) 

bCT k=l k k a 


The last equation yields the solution 


o = j/rN 


(43) 


with 


N 

J( C, N) £ 2 v,,* B * v, 

k=l 


T b ” 1 
k B 'k 


(44) 


where r is the number of measurements and N is the length of the data sample. The 

i a 2 

product rN is the degree of freedom in the quadratic form (44). The a estimate 
can now be used to scale the likelihood function; i. e. , 

N 

2 a 2 L = £ (v„ T B -1 v, + a 2 An det B + constants) (45) 

k=l k k 

In the F- 8 identifier, Equation (44) is evaluated at the minimum likelihood channel in 
order to estimate a 2 . This estimate is then used to scale all channel likelihood function! 
via Equation (45), before subsequent comparisons of these functions are made. v . 


Because the actual likelihood functions are being filtered, the degrees of freedom 
in the identifier are different from rN. They are determined by passing the constant 
r through the same filters and hence rise exponentially to the steady state level 
r/(l-u). 
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Likelihood Filters 


Figure 17 shows a high-pass filtering operation on the accumulated likelihoods in each 
channel of Algorithm 2. This operatiqnkeeps the accumulations current by exponentially 
de-weighting the pas fc, The rate at which de-weighting occurs (or the choice of time 
constants for the high pass) is determined by two conflicting requirements: 

1. A well-defined, correct minimum of the likelihood function requires slow 
de-weighting. 

2. Small tracking errors when aircraft parameters change requires fast de- 
weighting. 

The first requirement exists because likelihood functions do not necessarily have 
minima in the correct place (i. e. , at = c^) for short data samples. Some examples 
of - this behavior are given in Figure 19. The figure shows plots of experimental likeli- 
hood functions taken during our identifiability studies. The functions were accumulated 
without de-weighting and are shown at several time points. It is clear that times shorter 
than three to five seconds can give dramatically incorrect answers. It should be noted 
that this phenomenon is not inconsistent with maximum likelihood theory. The theory 
guarantees correct answers only asymptomatically. 

The second requirement is less esoteric. As the aircraft changes flight condition, 
its parameters change in ramp-like fashion. The accumulated likelihood functions 
will then be out of data and cause parameter estimates to lag behind the true parameter 
values. The faster we de-weight accumulated data, the less the lag. 

The design compromise for these requirements is a de- weighting time constant of five 
seconds. This satisfies the experimental curves in Figure 19 and produces 30 percent 
lag errors in M when the aircraft is accelerated with full afterburner starting at 
Flight Condition 5. (Traces for this maneuver are shown later in this section. ) Errors 
in M go of the order of 30 percent are quite tolerable for scheduling the control laws in 
Section 7. , It should be recognized, however, that the F-8C can be maneuvered so as to 
produce substantially faster flight transitions, and hence we must maintain tolerance 
in the control laws for errors well above 30 percent. 

In addition to low frequency de-weighting of accumulated likelihoods, it was also found 
desirable to prevent very high frequency data (such as sharp edged gust or command 
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Figure 19. Concluded, 



responses) from being accumulated. This was done by adding relatively high frequency, 
low-pass filter to the accumulation and high-pass network. The total likelihood filter 
then takes the following (digital) form: 



Unity gain low pass 

t = 0. 6 sec 

= exp (- A/0.6) 


Accumulation and 
high pass 

t = 5 sec 

p, = exp (-A/5) 


The symbols Ae k , Ae^, and e k denote generic inputs, intermediate states, and outputs, 
respectively. The filter is actually used to generate all likelihood functions, 
i = 1, 2, ... 5, and all components of VL and V L. 


CONTROL LAW ADJUSTMENT 

To complete the adaptive concept, the identifier described above was combined with a 
simple control gain calculation. Both pitch and lateral gains are adjusted as a function 
of the primary identified parameter, M gQ . In the pitch axis, the basic CCV control law 
gain schedule is utilized, but with a somewhat higher gain ratio. The CCV schedule 
(Section 7) calls for 

_ _ 0.158 

c* " — 

q 

This was replaced with the schedule 

_ 0.35 0.015 (47) 

vi . " — r — - ” a 


where the conversion from q to M go is based on the F-8C's very good linear q - M go 
relationship shown in the scatter plots of Appendix A (q = -23 M gQ ). The gain schedule 
(46) is a little more than twice as high as the CCV schedule, yet it still provides well 
above six db of stability margin. For additional protection, the schedule is restricted 
to a maximum gain output of 
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c* max 


0.35 

100 


0.0035 


and to a minimum gain of 


G 


c* min 


0.35 

‘600 


0.00058 


Similarly, the lateral gain was defined by the following schedule: 


g lat 



- 343 G c* 


( 48 ) 


This is automatically restricted to the range 


1.20 


12£ s 

100 ^LAT 


120 

600 


0.20 


by the limits imposed on G 

c ?l » 

SIGNIFICANCE TESTS 

Although the above schedules provide ample stability margin when the identifier is 
working properly, they do not exclude system instabilities for extreme identification 
‘errors such as might occur with no signals present in the control loop. To protect 
against this possibility, a self-checking feature was added to the identifier which monitors 
the statistical significance of M fio estimates. This self check calculates approximate 
theoretical probabilities that the true parameter values fall into the range covered by 
each channel. Whenever this probability exceeds a specified threshold (10“ 6 ) for any 
channel, the upper G ^ limit is reset to a gain value low enough to assure stability for 
true parameters actually in the range of that channel. The calculations associated with 
the self-check operation are summarized below. 

Probability Derivations 

Let Pr denote the conditional probability that the true parameters are actually in the 
range c/ 1 * covered by Channel i. This probability is defined by 
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( 49 ) 


Pr (l) = J n(1) P(C|Y n , D n ) d C 

p<y n Ic. yKii - 
V‘> P« N |U N >, 


Here p(C) is the a priori probability distribution for parameter vector c. This distribu 
tion can be assumed to be uniform ■with equal probability for each channel range; i. e. , 

p(6) = \ (i) for all C e 0 (l) 


X<i) J (i) « " ^ 


n v 


Then (48) reduces to 

(i) ' <i> 




p(Y N |U N ") p(Y N 1C< V ^(i) dC 


(51) 


P (Y N ic (i) . u N ) 
mp(y n |u n ) 


Next, it is convenient to divide through by Pr^ ^ “ 1 (i* denotes the min-L channel) 
and to take natural logs: 

r< i) 


Pr 


(i) o- ?£ 


(i) pffJr. V 


Pr®” p(Y N lt <i * ) .V 


(52) 


to Pr (i >-«nPr tt >-MPr (1 *> 


-L tt > + L«*> 
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This final equation forms the basis for the significance test. Since we want to detect 
(i) -6 

cases for which Pr' ' > 10 , it suffices tp test for 

L (l) - L (l * J :a -jen(10" 6 ) = 13.8 ( 5; 

and to reset G c ^ max whenever this test-fails , 

Gain Limit Reset 


The actual limit imposed on G in the event that the significance test fails is given by 


G 


c* max 


-0.039/1.5 M 



(54) 


This is eight db above the standard schedule. Equation (47), evaluated at the maximum 
I m 6q I Value covered by Channel i. Of course, if the significance test fails for several 
channels, then the smallest limit applies. 

DESIGN REFINEMENTS 


The identifier and control laws described above were tested on linear simulations in 
Minneapolis and on the nonlinear F-8 simulator at Langley Research Center with a 
high degree of success. They only posed one serious problem- -computing time. On 
the CDC-6600 computers at both simulation facilities, the identifier consumed approximately 
four to six msec per control cycle. Roughly scaled to the F-8C f s IBM AP-101 machine, 
this corresponds to 20 to 30 msec and exceeds the total frame time available at the nominal 
50 Hz update rate planned for F-8C flight tests. A serious effort was therefore under- 
taken to develop design refinements which would reduce running time of the identifier. 

Three such refinements were developed: 

1. Single-channel sensitivity calculations, 

2. Multi-rate structures, and 

3* Index-free Fortran code for primary channel calculations. 

These refinements reduce CDC computing times for the complete adaptive concept to 
approximately one msec/cycle. They are discussed in detail in the following pages. 
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Single Channel Sensitivities 


Since interpolation between channels of Algorithm 2 utilizes sensitivities from the min-L 
channel only, it is clear that a great deal of redundant computing is being done in the 
other channels. Redundant operations include sensitivity calculations for vx and likeli- 
hood accumulations for vL and v 2 L. These operations have been eliminated by providing 
a data hand-off procedure which allows one channel's sensitivity calculations to be 
started up and another's terminated whenever a channel change occurs. The transition 
occurs smoothly and without severe estimator transients. It is accomplished with the 
following transfer sequence from Channel j to Channel i at time t k : 

1. Initialize vx (t k > in Channel i 

Equations (39) with states q(0) and or(0) replaced by their current estimates 

a (i\ a (i) 

from Channel i, q (t k ), a (t k ). 

2. Transfer Accumulated v 2 L Matrix from Channel j to Channel i 

(55) 


, 1 , 1 } 


Matrix T adjusts the implied estimation accuracy of the transferred second - 
partials matrix. Basically, the accuracy is increased for every transfer in 
order to avoid rapid estimation transients. However, larger increases are 
used for transfers from high to low |M §q | values, and actual decreases are 
used early in a flight (through the parameter p). The larger increases preserve 
the identifier's constant percentage error property (Section 6), while the decrease 
speeds initial convergence. 
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V 2 L (i) = T 7 2 L (;i> T T 


where 


M 


(j) 


T = diag {p(l + 6 °/-r -) 

M ' ; 

6o 


with 


p = 1. 10 - exp (-t /5) 



( 56 ) 


3. Compute New VL Vector for. Channel i 


. (i) _ 2 t (i) . (i) - , 

VL = V L (c - £ ) 


This equation assures that- the identifier will interpolate to the current estimate 

A 

c k from the new channel immediately after transfer. Hence the parameter 
estimate remains continuous. 

After transfer, the propagation of Vx and accumulation of vL and V 2 L proceed as normal 
in Channel i, starting from the initial conditions V x (t^), VL^^, and v 2 L k ^ defined 
above. The same calculations in Channel 3 are terminated. 


Multi-rate Structure 


The second timesaving refinement of the adaptive concept consists of a rearranged 
sequence of computation. Since many operations do not need to be repeated at every 
control sample time, the calculations were broken into six separate groups, each 
containing those calculations which must be performed at every minor cycle plus a subset 
of calculations which can be performed less frequently. These six groups are then 
executed sequentially, one per minor cycle. The resulting multi-rate structure is 
illustrated in simplified flow diagram form in Figure 20. 

Fortran Implementation 

Other substantial time savings were achieved with modified Fortran code. In particular, 
a separate subroutine was prepared to implement five Kalman filter equations and one 
set of sensitivity equations in index-free Fortran. This subroutine accomplishes all 
minor cycle calculations shown in Figure 20, 

Figure 20 also illustrates several other features of the final identifier implementation. 
These include: 


1 * Temporary likelihood accumulation for each minor cycle in order to prevent 
data loss between likelihood filter updates (which occurs once per six minor 
cycles), 

a 2 

2. Limits on a to prevent excessively large noise level estimates induced by 
model mismatch at the min-L channel. 
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MLE Identifier Flow Diagram 
Figure 20* Continued. 









































3. Hysteresis on channel transfer decisions to prevent chattering when two channels 
have nearly equal likelihood functions, 

2 

4. Addition of small diagonal elements to V L in order to assure invertability, 
and 

A 

5. Limits on c to protect against gross estimation errors. 

MLE CONCEPT PERFORMANCE 

The complete adaptive control concept was verified with extensive nonlinear simulation 
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trials on NASA Langley*s F-8 simulator. System performance was assessed with 
respect to four "measures of goodness:" 

• Identification accuracy at fixed flight conditions, 

• Convergence characteristics, 

• Tracking characteristics for standard flight transitions, and 

• Response to major maneuvers and configuration changes. 

The controller was also "flown" under pilot control from the iron bird and under 
simulated relaxed- static-stability conditions. 

Accuracy at Fixed Flight Conditions 

Identification performance (and thus control gain accuracy) were verified by subjecting 
the closed-loop adaptive control system to a standard sequence of test conditions while 
in trimmed flight at various fixed flight conditions. The test conditions consisted of a 
period of quiet, followed by a 20 ft/sec 2 (6.1 m/sec 2 ) C*-command doublet, followed by 
a period of six ft/sec (1. 83 m/sec) atmospheric turbulence, followed by another command 
doublet. This entire sequence was run with and without sensor noise, and each run 
included a small random test signal inserted at the C*- command point. 

Accuracy tests were run at five flight conditions: FC1, FC5, FC8, FC10, and power 
! approach, FC17. Strip chart recordings for these flights are shown in Figures 45 
| through 54 in Appendix B. Each run is described by two pages of recordings. The 
first page shows major aircraft variables, and the second shows variables from the 
adaptive control law. 
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Performance summaries of the accuracy tests are given in Tables 15 through 18 for 
FC1, FC5, FC8, and FC10, respectively. These tables show estimation error percentages 
for the various test conditions. The errors were computed by comparing linearized 
derivatives (M gQ , M ff ) or actual quantities from the simulation (q, V, a) with estimated 
values. The estimates M^ # V, and q were taken directly from the identifier, 

and was synthesized as follows: 

A A 

(X , = a + 0. 5 
wL 

The summary tables show basic accuracies of 10 to 20 percent, which are very consistent 
with the predictions in Section 6. Poorest performance occurs at max-q "on-the-deck" 
(FC10), where errors and q errors are as high as 30 percent and velocity errors 
are as high as 40 percent. These errors are well within the stability margins of the 
control laws. 


TABLE 15. MLE IDENTIFICATION ACCURACY SUMMARY AT 
FIXED FLIGHT CONDITIONS 

FC1 




Parameter error (percent) 


Conditions 

q 

M 

6o 

M 

a 

V 

wL 

Quiet 

2 

0 

10 

4 

1 

Turbulence only 

8 

7 

2 

8 


Sensor noise only 

5 

4 

10 

15 


Turbulence plus 
sensor noise 

10 

7 

2 

1 


Pilot doublets 

2 

4 

2 

6 


Resolution of percent 
error calculation 

1 

4 

8 

3 

5 
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TABLE 16. MLE IDENTIFICATION ACCURACY SUMMARY AT. 
FIXED FLIGHT CONDITIONS 

FC5 



j 

Parameter error I 

percent) 


Conditions 

: q 

6o 

• M 

(X 

V 

°wL 

Quiet 

2 

5 

8 

8 

8 

Turbulence only 

5 

5 

8 

2 


Sensor noise only 

10 

5 

7 

17 


Turbulence plus 
sensor noise 

10 

5 

7 

3 


Pilot doublets 

4 

5 

7 

8 


Resolution of percent 
error calculation 

1 

ii 

15 

' 5 

5 


TABLE 17. MLE IDENTIFICATION ACCURACY SUMMARY AT 
FIXED FLIGHT CONDITIONS 

FC8 




Parameter error (percent) 


Conditions 

q 

M 6o 

M 

a 

V 

01 

wL 

Quiet 

16 

0 

2 

1 

4 

Turbulence only 

22 

6 

6 

1 


Sensor noise only 

15 

0 

2 

4 


Turbulence plus 
sensor noise 

30 

12 

2 

4 


Pilot doublets 

15 

0 

2 

4 


Resolution of percent 
error calculation 

1 

3 

2 

2 

16 





TABLE 18. MLE IDENTIFICATION ACCURACY SUMMARY AT 
FIXED FLIGfiT CONDITIONS 

FC10 




Parameter error (percent) 


Conditions 

n 

60 

M 

a 

V 

wL 

Quiet 

2 

2 

17 

40 

3 

Turbulence only 

20 

23 

9 

23 


Sensor noise only 

30 

29 

5 

29 


' 

Turbulence plus 
sensor noise 

15 

23 

9 

23 


Pilot doublets 

4 

6 

9 

29 


Resolution of percent 
error calculation 

1 

2 

4 

3 

9 


Convergence Characteristics 

2 

To verify convergence characteristics, the initial L, vL, and 7 L functions of the 
identifier were set to zero and the initial min-L index, i*, was set to a channel not 
covering the true flight condition* This was done for several combinations of channels 
and flight conditions and for various inputs* Strip chart recordings for these tests are 
given in Figures 55 and 56 in Appendix B, and the performance results are summarized 
in Table 19. Convergence times at or below one second were realized in all cases* It 
is particularly interesting to observe convergence characteristics obtained for C*-doublets 
without test signals (Figure 56, Runs 104 and 77). For these cases, the identifier 
drifts until the first command occurs and then responds almost immediately with proper 
estimates (and control gains). Comparing the initial transient with subsequent ones 
shows no major differences. This rapid response to commands may make it possible 
to operate the control loop without any test exitation. More experimentation will be 
required, however, to substantiate this potential. 
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TABLE 19. MLE CONVERGENCE SUMMARY 





Convergence 

Run 

Condition 

From 

channel 

To 

FC 

Time to 
80% 
(sec) 

71 

C*- doublets 

4 

5 

sSl 

68 

99% turbulence 

4 

5 

* 1 

103 

C*- doublets 

4 

1 

£ 1 

104 

C*~doublets 
no test signal 

4 

1 

a < 1 

77 

C*- doublets 
no test signal 

2 

1 

S < 1 

74 

99%' turbulence 

2 

1 

< 1 


a Measured from occurance of first command. 


Tracking 

To verify tracking capabilities, the aircraft was flown through two standard flight 
transitions: 

1. Max accelerations with full afterburner starting in trimmed flight at FC5, and 

2. Half-power deceleration starting in trimmed flight at FC8. 

Both transitions were flown with the CAS engaged (with zero C*-command). For the 
first transition, this produces rapid acceleration and a gentle climb to about 27,000 
to 28, 000 feet (8. 5K m) altitude and barely supersonic speed. The second transition 
produces slow decelerations to subsonic speeds with little change in altitude. 

Each transition was flown in quiet air, in 99 percent turbulence, and with 99 percent 
turbulence and sensor noise. All cases included a small C*-test signal. Strip charts 
for these cases are again given in Appendix B. Figures 57 through 59 cover the first 
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transition, and Figures 60 through 62 cover the second. Performance is summarized 
in Table 20 in terms of q tracking errors only (these were the easiest to calculate since 
q is available continuously from the simulation). 


TABLE 20. ^TRACKING ERROR SUMMARY 


Condition 

. 

Transition 1 

Transition 2 

Quiet 

- 31% peak 

+ 27% peak 

99% turbulence 

- 28% peak 

+ 38% peak 

99% turbulence 

- 28%, +32% peaks 

+ 38% peak 

plus sensor noise 


(+50% with speed 
brake applied) 


Maneuvers and Configuration Changes 


The adaptive controller was also tested under various maneuver conditions and configura 
tion changes. These runs are primarily intended to explore qualitative properties of the 
concept in large signal situations and to assure that no drastic upsets occur in normal 
operation. The runs are summarized below with figures given in Appendix B. 


Figure 62. 
Figure 63. 
Figure 64. 
Figure 65. 
Figure 66. 
Figure 67. 
Figure 68. 


Speed brake extended 

Roll maneuvers at FC1, FC5, FC8, FC10 
360 deg. rolls at FC1 
Step at gusts 

Gear-down transient at FC17 
Gear-up transient at FC17 
Wing-up transient at Mach 0.23, 3K 


No surprises were found in these tests. The identifier performs properly for standard 
roll maneuvers and alpha gusts. It developes tolerable 30 percent errors (in q) during 
360 degree rolls, 10 to 20 percent errors during gear transitions, and larger 50 to 70 
percent errors during wing transitions. None of these errors have major consequences 
on the controllability of the aircraft. In fact, when the wing transition was flown under 
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pilot control from the iron bird, the pilot judged control characteristics to be comparable * 
to the CCV laws. This same judgement was also rendered for the adaptive controller's 
performance over the rest of the flight envelope. 



The adaptive system was also flown with simulated aft motions of the c.g. For this 

purpose, the pitch control law was modified to include a low-passed pitch rate feedback 

2 

path developed specifically for RSS flight in the CCV program. This extra feedback 
is illustrated in Figure 21. It includes one additional adjustable gain, G R gg, scheduled 
as a function of M„ and M . The identifier was changed only in the nominal model ; 
parameterization. Equations (36), where the M # equation was replaced by M a = (0.3 +0.92 c 2 >* 
M 6q with larger a priori limits on c 2 (-± 0. 6). 

Flights were made at FC1 with zero and -2. 0 ft c.g. displacement. The latter is 
approximately neutrally stable. Transient traces are given Appendix B, Figure 69. 

The identifier is seen to have no difficulty with the large M a change induced by the 

c.g. displacement. ,, iV 



Figure 21. Relaxed-static-stability control law. 
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ML E PERFORMANCE ISSUES 

As shown by the nonlinear simulation tests above, the MLE adaptive control concept 
-performs well throughout the flight envelope* Its performance is not ideal, however, 
and it is important to highlight basic limitations. These include: 

• Tracking errors, 

• Effects of fixed statistics in the Kalman filters, and 

• Test signal requirements. 

The tracking limitation has already been discussed in the likelihood filter design section. 
We need a likelihood accumulation time of three to five seconds in order to assure 
correct minima of L, yet this accumulation produces lags when the parameters vary 
due to flight transitions or configuration changes. These lags are responsible for the 
largest errors seen in the performance trials. They can be relieved somewhat by 
additional cues to the identifier. For example, we could tell the identifier that the wing 
is going up or that velocity is increasing (forward acceleration) and allow it to adjust 
accordingly. Further research is recommended along these lines. 

Effects of the second limitation (fixed statistics in the Kalman filters) are evident in 
most of the accuracy runs. The estimates change somewhat between quiet or sensor- 
noise -only conditions and turbulence conditions. This is especially evident at high q 
(FC10), where the total change is as much as 50 percent. The problem is caused by 
filter gains which are designed for one condition (turbulence plus sensor noise) and 
then operated at another. It could be alleviated by providing additional data concerning 
operation conditions (i.e., turbulence estimates) to adjust the filter gains. Research 
is also recommended in this direction. 

Limitations due to test signal requirements are less quantifiable at this point. Test 
signal levels are primarily determined by pilot considerations and noise levels in the 
actual flight control system, so they will not be fully known until we fly the airplane. 

For the simulations at Langley Research Center, a small random test signal was found 
adequate. This signal used the spectral shaping network shown in Section 6 (with <« - 6, 

C = 1.25) and was inserted at the C*-command point. Its rms magnitude was adjusted 
to C* = 4.0 ft/sec (1.22 m/sec ), which produces the following approximate accelera- 
tions at the pilot station: 
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Test signal acceleration 

FC 

rms g*s 

m/sec 2 

1 

0.025 

0.25 

5 

0.010 

0.10 

8 

0.020 

0.20 

10 

' 

0.040 

0.40 


These signal levels can be lowered somewhat with tolerable (but noticeable) performance 
degradation. An example is shown in Figure 70 (Appendix B) which shows an accuracy 
test for FC1 with half the normal test signal. Potentially, they could be eliminated 
entirely if operation with rapid gain response to pilot commands proves acceptable. 

More experimentation is recommended to explore these possibilities. 
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SECTION 9 


THE MODEL TRACKER DESIGN 


This section presents the design of an adaptive gain adjustment procedure based on 
explicit parameter identification* The identified parameters are coefficients in a 
simple model of the linearized aircraft equations of motion. They are adjusted on-line 
to minimize a Liapunov error function. This section will present the theoretical frame- 
work of model trackers, discuss the design issues in a practical model tracker, and 
finally present the performance of a model tracker for the F-8C based on the pitching 
moment equation. The adaptive gain adjustment of the pitch and lateral CAS is made 
on the basis of an identified surface effectiveness parameter. 

PROBLEM FORMULATION 

The explicit model-tracking identification method is illustrated in Figure 22. The true 
system is represented in state variable form as 

. i ■ Fxt Gu (57) 



Figure 22. Model-tracking identification based on Liapunov function. 





A model is chosen to have the same structure 


x = Fx + Gu + v 
m 


(58a) 


where the elements of matrices F and G are adjustable. An error input to the model 

A A 

is denoted by v. Identification is complete when F = F, G = G, and v = 0. The procedure 

involves finding a way to adjust the parameters to force v to zero. 


A measure of the identification error is the vector difference between the system and 
the model states, e = x - x. A Liapunov design approach will provide expressions , 
for F, G, and v which force convergence and guarantee stability. * The input v 
is chosen to be v = D(x m - x) = De, where D is a matrix of design parameters. The 


resulting expression for the model tracker can be written as 


k * [F -D]x + Dx + Gu 
m m 


(58b) 


The error vector (e) defined above satisfies 
Is = De + (F - F)x + (G - G)u 


(59) 


The differential equations for F and G are derived using Liapunov's direct method to 
provide asymptotic convergence of e, F - F, and G - G to zero. One way to achieve 
this goal is to take the Liapunov function to be 


V ■ ~ [e'Pe + 2 p.zf] 
i . l l 


(60) 


with P as a positive definite symmetric matrix, {p.} a set of positive constants, and z^ 
components of the vector z of parameter differences; i.e., z' = (f^ - fjj, " f j2’ ** i ' m 
Then 


V = ~ e'(PD + D'P)e + S [p.Zj^ + h.(p, e, x, u)]^ 


(61) 


where h. is a sum of terms of the form e k Pkj(i) x q(i) or e k p kr(i) U s(i)* ChoOSing D SUCh 
that PD + D'P is negative definite and J - ' ,, 


z. = - p? 1 h. 
1 1 1 


(62) 


101 



yields a negative semidefinite ’O', which is zero only if e = 0. Then sufficient variation in 
u and x will yield convergence of the vectors e and z to zero. 


If P is a diagonal matrix, the h. simplify, and then the parameter adjustment given in 
(62) becomes a function of only one error and either a state or an input. Therefore, the 

A A 

individual elements of F and G should be adjusted using 


f., = X. .e.x. 
ij ij i J 


(63) 


g. . = u>..e.u. 
ij ij i ;j 


(64) 


where the "gains" and u>y must be positive numbers. 

Note that the Liapunov function could have been defined such that the e.'s in Equation (63) 

3 1 * . 

are replaced by odd functions of the e.'s such as sgn[e.], ore , etc. This will be 
considered later. 


Identification proceeds provided the system is being perturbed by some arbitrary input 
u(t) such that the parameter adjustments in Equation (63) are not identically zero. 

Recall the structure of Figure 22. By substituting the expression v = De, the form of 
the explicit model tracker shown in Figure 23 results. 



Figure 23. Model tracker structure. 
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DESIGN ISSUES 


Implementing the tracker of Figure 23 requires specifying parameters D, X^ and u^. 
Also, the states x and inputs u must be measured for use in the parameter adjustment 
equations. This section will present the design issues related to: 

• Modeling - - definition of the states and inputs in terms of measurable quantities, 

• Selection of D, X, and <u parameters, 

• Accommodating trim with the model tracker, 

• Effects of turbulence on the model tracker, and 

• Effects of measurement noise on the model tracker. 

Each of the above considerations influences the final design of the model tracker. 

Finally, the final functional block diagram and its F-8 simulator performance are 
presented. 

Modeling 

The F-8C aircraft models presented in Section 5 were simplified for the design of the 
model tracker. The linear model was reduced to the short-period perturbation equations. 
The effects of trim values will be discussed later. These can be written as 

q e M « + M q + M.6 (65) 

or q o e 

&= q + V + T a g +z 5 5 e <66> 

Fitch rate, q, and normal acceleration, N , are the measured variables where N is 

z z 

positive down and measured "d" m 'forward of the c.g. and is given by 

-N <* Z Vor + Z Vj. dfi (67) 

55 Or o 

(Here the notation follows the definitions of Section 6. ) 

To have the states be measured quantities, q and N are used as state variables. By 

z 

differentiating Equation (67) and using Equations (65) and (66), the short period can be 
written as 
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( 68 ) 


where the elements of F and G are obtained using a linear transformation as in Reference 2. 
Values of fjj„ fjg and g^ for the pitching moment equation are given in Table 21 for 
selected flight conditions spanning the F-8C flight envelope. 


TABLE 21. PITCHING MOMENT PARAMETERS 


FC 

h (ft) 
(Km) 

Mach 

c 

(psf) 

i 2 

(n/m ) 

f n 

f 12 

g ll 

1 

6.1 K 

0.67 

305 

14,603 

-0.7105 

0.3218 

-11.400 

5 

6.1 K 

0.40 

109 

5, 219 

-0.4446 

0.4817 

-5.024 

6 

6.1 K 

0.90 

551 

26, 382 

-0.9972 

0. 3017 

-21.980 

? . . 

12.2 K 

0.70 

134 

6, 416 

-0.3591 

0.3157 

-6.790 

8 

12.2 K 

1.20 . 

395 

.18,913 

-0.5466 

1.6020 

-19.070 

9 

3.0 K 

0.80 

652 

31,218 

-1 .1260 

0.2896 

-23.680 

10 

, 0 

0.70 

725 

34, 713 

-1.2650 

0.3200 

-25.230 

-ill-. 

0 

0.30 

133 

6, 368 

1 

-0.6204 

0.3566 

-6,032 

12 

0 

0.53 

416 

19,918 

-1.0060 

0.3405 

-15.860 

16 

12.2 

0.90 

222 

10,629 

-0.5704 

0.5802 

-11.040 

17 

F. A. 

0.18 

53 

i 

2,538 

. 

-0.3822 

0.1668 

-0.879 


The components of the input vector u (6 and 6 ) can be determined by including the 

e e 

differential equation describing the primary actuator from Section 5. 

: 6 * -12.5 6 + 12. 5 u (69) 

e e e 

The secondary actuator dynamics are assumed to be unity. Then u^ is the command 
to the secondary actuator. 

104 ; 



















Parameter 


Adjustment Equation 


The parameters D, X, and u> in the- model tracker were selected with the aid of a linear 
simulation of Equation (68). Since we wish the error to reach zero quickly, D was set 


equal to 




-5.0 




(70) ’ 


This is a time constant of 0. 2 second. 

The adjustment "gains" X.. and o>.. were picked to result in good parameter tracking | 

for step and ramp parameter changes. Values of X^ and optimized for a step | 

parameter change tended to be low for following a ramp parameter change. In both | 

■ ' ’ '• ' ■- ; * 
situations the values were dependent on the magnitude of the test signal (u e ). £ 

Starting values for the parameter gains were determined using small test signals. 

Step parameter changes and ramp parameter changes were made to gather a data base 
for further design work. Ramp parameter changes were generated by linearly inter- 
polating between FC5 and FC10. Although this does not reflect the actual way the | 
parameter changes, this approximation is sufficient for "tuning" the model tracker. » 

Typical parameter adjustment gains are shown in Table 22, The smaller values worked 
well for step parameter changes; the larger values were better for tracking parameter 
ramp changes. These values were initial starting points. As each design issue is 
considered, these values are modified, 

Some difficulty was experienced in achieving , "good" (monotonic) convergence of all | 
parameters simultaneously. This difficulty occurs even in the special case of tracking 
two elements in the first row of the system equations assuming all other parameters 
are known. This special case maybe written in the form 


*1 = f ll X l + f 12 X 2 + V<X 3* •••’ V U) 
y 1 = f u X l + f 12 x 2 + v(y 3 V u) + d(y i “ X 1 ) 

& j = Sf jk X k + Sg jkV i > 1 

where x is the true state and y is the model state. 


(71) 

(72) 

(73) 



TABLE 22. TYPICAL PARAMETER ADJUSTMENT GAINS 3, 



' *Vk + Es jkV J » 1 

e ■ y l ■ V *1 ■ V - hv 

The tracker eq uations take the form 

h = De + XjZj + x 2 z 2# d < 0 

. -1 •' 

Z l" " P 1 x l e * Pi > 0 

V - p’ 1 '*2®* p 2 > 0 


(74) 

(75) 

(76) 

(77) 

(78) 


This set of equations may be viewed as the following linear system with time varying 
coefficients: 
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with 



(80) 


The matrix A has two eigenvalues with negative real parts and a zero eigenvalue for 
each value of x 1 and x 2 . Thus, it is not surprising that for specific x^ and x g time 
histories, the convergence of and is initial condition dependent. In some cases 
one parameter will converge essentially monotonically as the other exhibits initial 
divergence. To alleviate this difficulty, the following nonlinear tracking equations were 
studied: 



(81) 


where f° and f° 12 denote nominal (constant) values of f^ and f 12 , respectively. The 
convergence of this system is assured via the Liapunov function 

V 8 e 2 + pjZ J + P 2 z| (8D 

if p x >0, p 2 > 0, « 0, and D< “ f ° x > 2 + P 2 ^2^12 “ f 12 )2 ^ 4 * 

A priori bounds may be used for (f^ - and (f 12 - f° 2 ) 2 in the constraining inequality 
on D. 

After some experimentation, it appeared that only limited improvement was noted, 
and even this required f?. close to f... As a result, no changes were made to the 

i 3 13 

parameter adjustment equations. 
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Trim Considerations 


As explained earlier, only perturbation equations have been considered. The trim 

values of the accelerometer and elevator surface will now be included. An additional 

parameter (M^) representing the trim pitching moment was added to the model where 

N and 6 now have non-zero trim values,* .The adjustment equation used for estimating 
z e 

M was 
o 

M = - X e (82) 

o 

A 

Simulation runs were made to pick a X such that M q would account for changing 
flight conditions. Early results show that tracking M q degraded the tracking of the 
remaining parameters. The M q parameter was dropped, and high-pass filters were 
added. to remove the trim values of N and 6 to permit tracking the perturbation equations. 

Z 6 

Early results with first-order high-pass filters were favorable. However, in tracking 
ramp parameter changes, large biases appeared as a result of the high pass "differentiating" 
the slowly varying elevator trim. This was alleviated by using second-order high passes. 
There is a trade-off in selecting the high-pass time constants. The time constants were 
chosen to compromise rapid trim changes (such as roll maneuvers) with good sensitivity. 

A one-second time constant is currently used. Finally, the same filter is applied to the 
pitch rate gyro even though its trim value is zero. This is done to maintain the proper 
phase relationship between q, N , and 6 . 

Effects of Turbulence 


The effects of «_ (refer to Equation (68)) driving the true plant (but not measurable by 
the model) were determined using a linear simulation of the short-period dynamics. 

The angle-of-a ttack gusts were modeled using a first-order approximation to the Dryden 
spectrum as explained in Section 5. 

A simple C* feedback controller (with fixed gains) defined u g (t), and no additional 
pilot inputs or extra test signals were used. 

Simulation results indicated that identification of the first row of the model (f^. 
f 12 , and g n in Equation (68)) proceeded satisfactorily, but the second row did not 
appear to converge or converged to wrong values. 
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This is not unexpected since the model does not contain a gust input. If angle -of -attack 
were measured, then this problem would be alleviated. However, since we do not 
wish to use air-data measurements (sucn as or)/ it does not appear that good tracking 
during gust disturbances can be obtained for the Coefficients of the normal force equation 
(the second row). An alternative would be to try to estimate the wind gust in the model. 
However, the identifiability results indicate that this cannot be estimated with sufficient 
accuracy. A possibility for tracking second-row elements would be to use thresholds 
to eliminate gust effects and to attempt to track only during large maneuvers. This ; 
was considered too restrictive for further consideration. j 

• ■ ■ ■ • • . ■ ' • " ‘ i : • 

Measurement Noise Effects ! ; 

Sensor noise was added to the q, N z » and 6 g measurements, and errors in tracking the 
first-row parameters (fj j, f 12 , and g^) were studied. The sensor noise was defined 
in Section 5. Initially, sensor noise seriously degraded parameter estimation as usually 

fiQ . . 

results in Liapunov designs. This is not surprising since the above errors are of the 
same order of magnitude as the vehicle motion induced by gust • or small test signals'. 

Several modifications to the tracker were made to reduce errors caused by sensor noise. 
This included filtering the measurements and adding thresholds and limits to the error 
function. 

First-order low-pass filters were added to remove as much of the high frequency energy 
as possible without degrading performance in the noise-free case. The time constant 
of the low pass was set equal to the time constant of the elevator primary actuator. 

Next, threshold and amplitude limits were added to the error quantity to reduce the : V 
amount of parameter adjustment resulting from sensor noise. The nonlinear element 
shown in Figure 24 was added to the error equation. The parameters Tj and T g were 
determined from simulation results. Tj is the error threshold and was set to 0. 5 x 10 
rad/sec. Mismatch in pitch rate less than Tj will not cause any parameter adjustment. 
The error is limited to T 0 = 0.0055 rad/sec to prevent large commands and disturbances 
from producing large parameter adjustments. 

The parameter gains that work well were 
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* 100 

V 5 ' (83) 

«j » 2. X 10® 

Further reduction in sensitivity to error magnitudes (especially for thunderstorms 

and large commands) was achieved by adding a second threshold element to alter X 

and «>. The parameter adjustment gains are reduced by a factor of 20 when the bandpassed 

elevator servo signal exceeds 0. 002 radians. 

Finally, rate limits were added to the three parameter adjustment equations to further 
equalize performance with small and large signal levels. They were chosen to provide 
adequate tracking rates for realistic maneuvers, but not to permit excessive parameter 


adjustment activity. The rate limits are 


i*u[ 

£ 

4„0 


I'lll 

£ 

0„1 

(84) 

ll n l 

£ 

0„1 
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GAIN ADJUSTMENT 


A 

The G c * gain used in the pitch CAS is determined from g^. A low-pass filter was 
added to smooth the estimate, and limits were added to constrain the gain values during 

A 

large transients in g^. Figure 25 shows the variation of gjj with flight condition 
plotted as a function of dynamic pressure. When the corrections for flexibility are 
included, an approximate fit is 

-g u = 0.0378 q + 1.2 < 85 > 



100 200 300 400 500 600 700 

(4.8) (9.6) (14.4) (19.2) (23.9) (28.7) (33.5) 


q, psf (kn/m2) 

Figure 25. Parameter g n versus dynamic pressure. 


Ill 



(86) 


Using the estimate of g^, dynamic pressure was estimated from 
q = 26.455 (-g^) - 31.75 

The error in dynamic pressure 

A 

kzJL 

q 


(87) 


is one of the time history plots included in the simulation data and used for accuracy 
evaluation only. 

FUNCTIONAL BLOCK DIAGRAM 


The final form of the pitch axis model tracker is shown in Figure 26. It contains a 

perturbation pitching moment equation with three parameters adjusted on the basis of 

a Liapunov error function. The filtering on q, N , and 6 is shown. 

z e 

The thresholds, nonlinear error functions, and rate limits previously discussed are 
defined. The specific values of these nonlinear elements were chosen on the basis of 
simulation responses. Finally, since the servo position rather than actual surface 
deflection is measured, a first-order linear model of the primary actuator has been 
incorporated. 5 

DIGITAL MECHANIZATION 

The tracker was designed using the continuous form of the equations and a Runge-Kutta 
integration method. 


The tracker was converted to a 32 sps digital version for incorporation with the digital 
control laws presented in Section 6. The filters were converted to a discrete equivalent 

• A A A 


using Tustin's method. The q, f^, 
using an Adams -Moulton formula. 


f 12 , and gj nonlinear equations are integrated 


X k + l = \ + l (10 \- 5 \-l + V2 ) 

with A = 1/32 second. 
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A flow chart of the digital model tracker is shown in Figure 27. 

PITCH AXIS PERFORMANCE 

The model tracker was evaluated against three "measures of goodness: 11 

1. Accuracy at several fixed flight conditions including effects of servo: 
hysteresis, sensor noise, turbulence, and pilot commands, 

2. Convergence properties, and 

3. Tracking varying flight conditions. 

The major performance measure is the accuracy of identifying the surface effective- 
ness parameter. Time histories from the NASA/LRC F-8 simulator for the conditions 
summarized below are contained in Appendix C. 

Accuracy at Fixed Flight Conditions 

The results of the accuracy runs are tabulated in Table 23. A standard pitch sequence 
consisting of periods of quiet, gusts, and C*-doublet commands was used as system 
inputs. 

The accuracy in identifying 'g^ varies from five to 40 percent for pilot command and 
turbulence. Sensor noise (only) causes errors to grow as large as 80 percent. 

Convergence Properties 

Convergence characteristics of parameter estimates are important in the model tracker. 
This test determines how well initial condition errors are removed when turbulence or 
pilot commands are applied. Two cases were studied. First, the algorithm was 
initialized at FC10 with the aircraft at FC5. Second, the algorithm was initialized 
at FC5 with the aircraft at FC1. Two separate inputs were applied for each case. The 
first input was turbulence of the Dryden form. The second input was C* pilot commands 
Results are summarized in Table 24. 










TABLE 23. PERFORMANCE AT FIXED FLIGHT CONDITIONS 


Input condition 

'• A 

Error in g^ (percent) 
Flight condition 

Comments 

■ 1 

5 

8 

10 

Quiet 

0 

0 

0 

0 

During quiet periods, no parameter adjustment is 
made; so if initial conditions are correct, no 
errors will occur for no test signal case. 

Pilot doublet 

5 

20 

30 

5 

- 

Turbulence 

10 

20 

40 

5 

- 

Doublet plus turbulence 

15-20 

10 

10 

5 

- 

Sensor noise 

20 

70-80 

25 

5 

Sensor noise causes parameters to drift off. 

Doublet plus sensor 
noise 

5 

10 

5 

5 

Doublet provides information that reduces errors 
caused by sensor noise only. 

Turbulence plus sensor 
noise 

20-25 

20 

40 

5 

Sensor noise effects dominate. 

Doublet plus turbulence 
plus sensor noise 

10 

10-20 

10 

5 

- 

Servo hysteresis 

70 

- 

i 

- 

- 































TABLE 24. CONVERGENCE RESULTS 


Test condition 

Gusts 

Commands 

Case 1 
IC = FC10 
Plant = FC5 

A 

gf l to 50% of steady state in 2 sec 

A 

g n to 80% of steady state in 4 sec 

Same as 
gusts 

Case 2 
IC = FC5 
Plant = FC1 

A 

gn to 80% of steady state in 4 sec 

Same as 
gusts 


Parameter Tracking 

Two profiles for changing flight conditions were run. The first started at FC5 (q = 5.03 

kn/ m 2) and accelerated at full power holding a constant value of C* (no pilot input). The 
aircraft accelerates from 61 m/sec to 152 m/sec and climbs 1524 m in about a minute. 
This profile was repeated for a quiet case, turbulence, and turbulence plus sensor noise. 


The second profile consisted of a deceleration from FC10. Thrust is reduced to idle, 
and the speed brake is extended. Altitude remains nearly constant* and airspeed drops 
off. Results are summarized in Table 25. Because of the thresholds in the tracker, 
some error buildup is necessary before any parameters can change. 


LATERAL AXIS PERFORMANCE 

The lateral- directional CAS used the structure described in Section 7 with the loop 
gain adjustment slaved to the pitch axis tracker. The gain G LAT I s computed from 
estimated dynamic pressure. Equation (86). Performance of the lateral-directional 
CAS in terms of turn coordination and Dutch roll damping is similar to the performance 
of the MLE lateral-directional CAS, so time history traces will not be repeated in this 
section, Performance also compares with the results presented for the high-gain 
limit-cycle design in the next section. 






















SECTION 10 


HIGH-GAIN LIMIT-CYCLE DESIGN 


This section discusses the design and performance of the limit-cycle gain changer (LCGC) 
applied to the F-8C. The relationship of LCGC to the aircraft CAS is shown in Figure 28. 
The gain changer is a closed-loop controller that adjusts the loop gain of the CAS based 
on secondary servo activity. The objective of the gain changer is to keep the CAS 
operating at the highest practical gain. In the absence of gusts or measurement noise, 
it will maintain a small amplitude limit cycle on the secondary servo. Adjustment 
of the CAS loop gain by this method is, then, an implicit adaptive process. This concept 
is similar to a model reference implicit concept since the objective of both approaches 
is to have the bandwidth of the aircraft's control system sufficiently high so that it can 
follow commands generated by an explicit model. In the case of the F-8G, the pitch 
axis uses an explicit second-order C*-model, and the roll axis uses an explicit first- 
order model of roll rate. 


The LCGC concept is more restrictive than implicit model reference approaches since 
the structure is limited to adjustment of a single gain parameter (i. e. , "loop gain' ). 
However, for the F-8C, a single loop gain provides sufficient adaptation. 



Figure 28. LCGC operation. 
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EVOLUTION OF THE LIMIT-CYCLE ADAPTIVE CONCEPT 


The limit-cycle concept evolved from early work (1957-1958) in adaptive controls 
sponsored by the Wright Air Development Center (Reference 70)* This work studied 
application of a bi-stable element (relay) to high -bandwidth feedback control loops. 

The results were generally unsatisfactory because acceptable limit-cycle amplitudes 
could not be achieved with adequate control authority. The concept of varying a linear 
gain to maintain the limit cycle at a fixed small amplitude was then conceived, thereby 
avoiding a constraint on control authority. The limit cycle was sensed via a bandpass 
filter, compared to a set magnitude reference, and the resulting error applied to a 
proportional-plus-integral function to vary the control-loop gain. The first application 
of this approach was an experimental system in an F-101A airplane (1958-1959), involving 
adaptive gain control in pitch and roll control augmentation systems (Reference 15). 

This was followed by the design, build, and flight test of a three-axis augmentation 
system for one of the X-15 aircraft (References 8 through 14) in 1964. The X-15 applica- 
tion used the adaptive gain control to maintain high -bandwidth loops for the aerodynamic , 
controls and also used the resulting gain level to signal engagement of reaction controls 
when the normal surfaces became ineffective. 

Both the F-101 and X-15 applications encountered higher frequency structural modes 
which reduced attainable gains below initial design levels. Notch filters were added 
'in each case; to adequately restore loop bandwidths. Some reduction in limit-cycle 
frequencies occurred as a result of the added filtering, necessitating some changes in 
parameters of the adaptive gain controls (e.g. , bandpass filter frequencies). 

Another problem common to the adaptive approach described above is the effect of 
system disturbances with significant signal content around the limit-cycle frequency 
(e.g., turbulence, sensor noise, high frequency pilot inputs). 

The gain control can interpret these as excessive limit-cycle amplitude and drive the 
loop gain down. Although generally satisfactory operation can be achieved for a stable 
vehicle, an unstable vehicle will generate added oscillations as the gain reduces, causing 
complete loss of control. To correct this problem and to improve gain regulation for 
a stable aircraft, a low frequency signal path phased to cause a gain increase was added 
to the adaptive controller. This new approach was never flight tested but was extensively 
studied, and hardware was built for launch boosters (Reference 71) and other aircraft 
(Reference 16). The design for the F-8C includes this "up-logic" feature. 
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A significant benefit of the above adaptive developments was experience gained with high- 
bandwidth control systems. By accommodating much of the airframe variations through 
normal feedback control, avoidance of any control parameter changes is sometimes 
possible or requires, at most, use of simple air-data schedules. Ironically, this 

situation eliminates some of the original motivation for adaptive development. 

GAIN CHANGER OPERATION 

The basic elements of the gain changer as applied to the pitch and lateral-directional 
CAS are shown in Figure 29. In operation the set point will increase the loop gain until 
the secondary actuator begins to oscillate at the limit-cycle (or 1 crossover ) frequency. 

The down-logic portion of the gain changer senses the resultant servo motion and decreases 
the gain until a stable limit cycle is established at the "linear critical gain. " Some 
proportional gain is required to stabilize the gain changer loop; otherwise it will tend to 
oscillate about the critical gain. However, if the proportional gain is too high, excessive 
gain changer activity results. Hence a compromise must be reached. In a quiescent 
environment, da m pin g of the loop gain response resulting in an overshoot of 100 percent 
in the servo limit-cycle envelope is generally acceptable. In actual operation, usual 
| noise levels will essentially eliminate this overshoot. The limit-cycle amplitude is 
j determined by the ratio of set point to down-logic gain. The down-logic gain is adjusted 
to give the desired limit-cycle servo amplitude (i0. 0035 rad peak-to-peak was used 
for the F-8C). Subsequent adjustment of limit-cycle amplitude may be made with the 
set point parameter. - 

DESIGN ISSUES 

The major issues involved in the design of a limit-cycle gain changer concern: 

• Selection of limit- cycle frequency, 

• Selection of limit-cycle amplitude, 

• Design of down-logic and up-logic bandpass filters, 

• Requirements for nonlinear LCGC response, 

• Definition of the set point, 

• Rate limits on gain adjustment, and 

• Amplitude limits on gain adjustment. 
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Figure 29. Basic elements of limit-cycle gain controller (LCGC). 









Each of these topics will be individually discussed. Next the functional block diagram of 
the LCGC as evaluated on NASA/LRC*s F-8C simulator will be presented. Finally a 
summary of the performance of the pitch axis LCGC and the roll axis LCGC will complete 
this section. 


Selection of Limit-Cycle Frequency 

The limit-cycle frequency is selected primarily to reflect the vehicle property being 
compensated for by the adaptive gain adjustment. Other secondary criteria include 
minimum interaction with normal control dynamics (a limit- cycle frequency higher than 
dominant frequencies), small variations due to flight condition and component tolerances, 
and acceptable amplitudes (avoidance of low-damped resonances). 


The main object of the gain adjustment applied to high-bandwidth aircraft control loops 
is to compensate for varying surface effectiveness (fuselage rotational acceleration 
per unit control deflection). Measurement of this quality with a single frequency signal 
must be performed at frequencies above aerodynamic roots (e.g., short period, dutch 
roll, or roll subsidence), but below structural modes (flexure). Fortunately, this also 
generally satisfies the secondary criteria listed above. With a limit-cycle frequency 
so selected (generally around three to four Hz in a fighter- type airplane), significant 
system dynamics usually include the rigid aircraft (a pure inertia with angular rate 
feedback), the power actuator (a first-order lag), the secondary servo (a second-order 
lag), and applicable system filtering (perhaps a notch for flexure mode attenuation). 
These dynamics are relatively stable, resulting in a predictable limit-cycle frequency. 
Applied to a pitch rate loop with no added filtering, the applicable transfer function is 




(—=— 
T S+l 
a 


) < 


■ + 2£a)S + 



where 

q = pitch angular rate. 

Mg = pitch surface effectiveness, 

e 



t = actuator time constant (e. g. # 0.05 sec), 
a 

go) = servo damping and frequency (e.g., 0.7 and 60 rad/sec), and 

6 - elevator deflection, 

e 
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Note for this example, a limit-cycle frequency of about 27 rad/sec would occur* The 

critical control-loop gain (K ) would vary inversely with M , approximated by 

c 6 e 

K M 

c 8 e . 

1 1 

T a (27) 2 

Limit Cycle Amplitude Considerations 

The amplitude of the servo limit cycle is designed to be twice the assumed value of flow 
hysteresis in the secondary actuator* Depending on the actual level of hysteresis 
or other control system nonlinearities , it may be necessary to adjust the limit- 
cycle amplitude with the set point parameter. The assumed amount of hysteresis will 
reduce the limit-cycle frequency by about 10 percent, which will change the limit-cycle 
amplitude. Final adjustment with actual hardware may be required* 

For the F-8C, the nominal limit-cycle amplitude is set to provide a peak-to-peak 
secondary actuator amplitude of 0.0035 rad. This is equivalent to 0.0021 rad elevator 

position and 0„ 17 degrees aileron position with the assumed power actuator characteristics* 

_ 2 2 
At FC1 (q - 14603 n/m ), this size of limit cycle produces o^ily 0.047 nf m (rms) normal 

acceleration. 

Down-logic and Up-logic Frequency Characteristics 

For the F-8C, the roll axis crossover frequency was intentionally set to be identical to the 
pitch axis crossover frequency at a value near 20 rad/sec (three Hz). The bandpass 
filters for the pitch axis are shown in Figure 30 and for the roll axis in Figure 31. 

There are only minor differences in the gain value. 

Rapid attenuation outside the frequencies of interest is desirable to minimize the effects 
of wide spectrum noise or disturbances on gain changer operation. Disturbances usually 
result in a reduction of loop gain which may degrade the performance of the CAS. The 
bandpass networks were implemented as third-order linear filters, resulting in a 60 db/ 
decade attenuation outside of the bandpass frequency. 

The down-logic bandpass filter is peaked slightly above the crossover frequency in each 
axis. This assures limit cycling on the positive slope of the bandpass, where phase 
lead will contribute to stability of the limit cycle and give more positive control of the 
limit -cycle amplitude. 
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The up-logic counters the tendency of the loop gain to reduce during gusts and other 
disturbances including pilot commands. The up-logic bandpass and gains maintain an 
energy balance during wide spectrum disturbances. Based on past experience, the up- 
logic frequency in this application was set a decade lower than the down-logic. If a low 
frequency instability exists (as in some previous limit-cycle applications), then this 
frequency would govern. Up-logic gain is set to maintain an average six-db gain margin 
during turbulence or pilot command inputs. Further optimization of the limit-cycle 
concept will require close evaluation of the up-logic mechanization, particularly during 
pilot inputs. Availability of spectral density information on pilot flying tasks would 
allow more accurate analytical definition of up- logic parameters. 

Nonlinear Gain Element 

The purpose of the nonlinear element is to maintain constant gain changer dynamics 
at all flight conditions. Near critical gain, a given percentage change in loop gain results 
in a proportional rate of change in limit- cycle amplitude at any flight condition. Therefore, 
*the resulting feedback (i.e,, the sensed limit-cycle amplitude) to the gain changer should 
produce a fixed percentage change in control gain. Since an exponential gain element 
has a slope proportional to the output (which is the control gain level), the above require- 
ment for a nonlinear element is satisfied with an exponential function. For the F-8C, 
the exponential function has a 50-to-l dynamic range and was scaled to produce a 
maximum value of unity. A 50-to-l range in loop gain will allow the gain changer 
to maintain a limit cycle throughout the flight envelope, except power approach. At 
power approach, the maximum value gain is insufficient to sustain a limit cycle. The 
maximum value of is limited to 0. 01. 

Set Point Definition 

Starting from a stable gain, the initial rate of gain adaptation is determined by the set- 
point and the integral gain. The set point and the integral gain are set such that their 
product will result in the driving of the nonlinear element from maximum to minimum 
(and vice versa) in a length of time consistent with the aircraft flight condition change 
capability. In this F-8C application, 40 seconds was selected as the nominal value. 

From other considerations, the set point and integral gain may require redistribution, 
but their product has been defined. 
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The integral rate may be further optimized when the system is subjected to disturbances 
to avoid excessive gain activity, particularly when the error limits in the system are 
reached. 

Rate Limits on Gain Adjustment 

Limits are imposed within the gain changer in order to define and bound gain changer 
operation when inputs are sufficient to saturate various elements. An important 
consideration is gain changer response to pilot inputs. A system operating at critical 
gain will obviously oscillate when disturbed. Ideally, a system with six-db gain margin 
provides the highest practical gain without excessive control surface oscillations 
(ringing). The rate limit should be set such that the product of the rate limit and 
proportional gain will produce a six-db change in the output of the nonlinear element. 
Thus, a rapid gain reduction will occur for pilot or disturbance inputs. Selection of 
the rate limit must also consider initial stabilization of a statically unstable vehicle 
if initialized at an unstable gain. 

A seoond vital consideration in specifying the rate limit is related to pilot activity 
at the down-logic bandpass frequency. If the pilot concentrates an input at this frequency, 
the gain could be driven down from maximum to minimum in about 12 seconds if the 
rate limit is saturated. This level of determined, sustained effort is not felt to be 
probable but should receive investigation in pilot-simulation studies effecting not only 
the rate limits, but the up-logic bandpass and gain. 

Amplitude Limits on Gain Adjustment 

Down-logic is given twice the authority of the up -logic in the adjustment of the auto- 
pilot gain (refer to Figure 32). The down- logic authority limit is set equal to the rate 
limit. For saturation in the up- or down-logic, the rate limits will determine gain 
adjustment. With both up- and down-logic saturated, the relative authorities assure 
that the gain will be reduced. Without this assurance, extreme noise environment 
could conceivably drive the gain above critical values. In the present studies, no 
single disturbance could create this situation; but possibly multiple disturbances, 
including pilot inputs, may although the possibility is considered remote. 



FUNCTIONAL BLOCK DIAGRAMS 

The final design of the pitch axis high-gain adaptive controller is shown in Figure 32* 
The majority of the parameter trade-offs in the design of the high-gain limit cycle 
were made using a two-degree-of-freedom analog simulation. The digital implementa- 
tion of LCGC used the Tustin equivalent of the bandpass filters and signal shaping. 

For this concept, the C*- model was changed slightly by adding a 0.1 second first-; 1 
order lag to the pitch stick shaping to reduce elevator "ringing" in response to stick 
commands. 

; in verifying the design on the Langley nonlinear simulation, it was necessary to add 
small lag <t = 0. 02 sec) ahead of the elevator servo to obtain the desired pitchaxis 
limit- cycle frequency. It is felt that this extra phase lag compensates for approxima- 
tions in the digital integration used in Langley's simulation. It probably would not be 

needed in the real hardware system and is not shown in Figure 32, , 

The roll axis of the high-gain limit-cycle adaptive system is shown in Figure 33. The 

gnin G ,, is set at 5.0*G T A _ and the roll axis is held near critical gain by adjusting 
6 roll LiA 1 

°LAT *** 016 lateral LCGC * The y aw a ^ ±s s ain “sed g lat ldlr ® cfl y* and hence thi& 
adjustment is slaved to the roll axis. The G LAT gain range for the yaw controller 
provides good performance and sufficient stability margins over the flight envelope. 

PITCH AXIS PERFORMANCE 

The high-gain limit-cycle design was evaluated against three areas related to adaptive 
performance: ;; ‘ r \ 

1. Accuracy at fixed flight condition (consider effects of servo hysteresis, 
sensor noise, turbulence, and pilot commands), 

2* Convergence properties, and 

3, Tracking varying flight conditions. 

Each of these areas will now be discussed. . 
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Figure 32. F8-C high-gain limit-cycle pitch axis . 





Roll Stick 



(Up) 


Figure 33. F-8C high-gain limit-cycle lateral- directional axes. 




Accuracy at Fixed Flight Conditions 

The standard pitch sequence consisting of periods of quiet, turbulence, and C*-doublet 
commands were used to measure the performance of LCGC and Flight Conditions 
1, 5, 8, 10, and 17. 

The time histories are contained in Appendix D. They were obtained from the Langley 
’ nonlinear simulation and were consistent with traces obtained on a linear two-degree - 
of-freedom analog simulation at Honeywell. The accuracy results are summarized 
in Table 26. Overall summary comments are provided together with the gain values 
(expressed as percent of critical) as functions of input conditions. 

Convergence Properties 

The concept of Converging" to a gain value does not apply to the LCGC in the same 
way it does to the explicit identifiers. In normal operation the gain changer will not 
be at a gain higher than critical. If the system starts at a value below critical, the 
rate of gain increase is largely determined by the integral gain. For example, if 
the LCGC is initialized at FC10 with the aircraft at FC5, it takes on the order of 10 
seconds (in gusts) to reach "steady state" (which is six db below critical). Pilot 
commands do not speed up "convergence" since their effect tends to decrease the gain. 

Adaptation to Changing Flight Conditions 

In the case of increasing and decreasing M$ , the gain changer tracked 
gain (quiet case) or tracked the "turbulence" gain value. 

Deceleration profiles starting from FC10 by reducing thrust to idle and 
speed brake show that the limit-cycle amplitude remains constant. 

Sens or N ois e Eff e c ts 

The LCGC was evaluated with the measurement noise defined in Section 5 for the 
pitch rate gyro, normal accelerometer, and servo position transducer. While 
gain changer activity increased in the presence of sensor noise, only the pitch rate 
gyro noise affected the LCGC to any significant degree. The gain changer will attempt 


the critical 


extending the 
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TABLE 26 i. PERFORMANCE AT FIXED FLIGHT CONDITIONS 



Gain (percent of critical) 




Flight condition 



Input condition 

1 

5 

8 

10 


Comments 

Quiet 

100 

100 

100 

100 

(1) 

Typically, the set point increases the gain 
over its critical value until the servo 
oscillation builds up. Then the gain is 
reduced to the critical value with a limit - 
cycle overshoot of about 100 percent. 

Pilot doublet 

60-63 

60-75 

60-71 

46-48 

(2) 

Gain is self-adjusting to achieve a nominal 
six -db gain margin during pilot inputs. Re- 
sponses have typical six-db "ring. " Up-logic 
prevents successive inputs from driving gain 
to minimum, establishing an average gain 
higher than six db. 

Turbulence 

75 

56 

81 

69 

(3) 

Gain is self-adjusting to provide control 
surface energy, induced by turbulence, to be 
essentially equivalent to the nominal limit- 
cycle energy (slightly higher due to up- 
logic bandpass). The net effect is a desirable 
ride softening as a function of turbulence. 

Doublet plus turbulence 

75 

56 

81 | 

69 

(4) 

Pilot inputs during noisy operation will cause 
gain to establish an average value about equal 
to sensor noise only case. Pilot input will 
initially drive gain down, wherever the 
starting point is. 

Sensor noise 

30 

20 

43 

48 

(5) 

Same as Note 3 except gyro noise reduced 
gain more than turbulence. 

Doublet plus sensor noise 

30 

20 

43 

48 

(6) 

See Note 5. 

Turbulence plus sensor 
noise 

30 

20 

43 

48 

(7) 

See Note 3. 

Doublet plus turbulence 
plus sensor noise 

30 

20 - 

43 

48 

( 8 ) 

See above Notes 3, 5. 

Servo hysteresis 

74 

68 

67 

80 

( 9 ) 

Frequency reduced 10 percent. 









to adjust the gain downward until the noise-induced servo motion matches the equivalent 
limit-cycle energy in the down-logic network* The noise sensitivity was studied by 
simulating three different levels and bandwidths of gyro noise (white noise, first- 
order filter). 

A. 0,0087 rad/sec rms with a two Hz bandwidth. (This value was suggested by NASA 
from flight test recording data.) 

B. 0,026 rad/sec rms with a 16 Hz bandwidth. (This is a conservative value defined 
for these studies.) 

C. 0.0035 rad/sec peak-to-peak with a 16 Hz bandwidth. (Obtained from measurement 
of actual Honeywell GNAT gyro is a pristine environment.) 

It is felt that actual gyro characteristics will fall somewhere between B and C above, 
in which case, the adaptive gain will tend to be between six and 12 db below critical. 
Significantly, this gain level is still greater than the original CAS q-schedule so auto- 
pilot performance should remain at acceptable levels. 

Turbulence Effects 


Thunderstorm turbulence will produce about the same gain reduction as gyro noise. 

At the 99 percent turbulence level (exceeded four days /year), the gain reduction is 
about three to six db, the greater reduction occurring at low q. Up-logic plays a role 
in maintaining the gain level. Actually, the gain reduction, especially for high turbulence, 
produces a favorable "ride softening" effect. 

Effect of Pilot Inputs 

Since previous limit-cycle systems have been prone to excessive gain reduction/^ 
special note is made of gain levels during periodic pilot inputs. Initial response to a 
pilot input is a six-db gain reduction, a desirable and necessary feature of gain changer 
operation. However, without some up-logic capability, successive pilot inputs tend to 
further reduce the gain. The resultant interplay of the up- and down-logic results in 
considerable gain activity, but the average gain is generally slightly less than 50 percent 
below critical. When the up-logic was removed as an experiment during the pilot 
input sequence, the gain quickly reduced to 20 percent of critical. 

f 
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The effects of sensor noise, turbulence, and hysteresis on the autopilot gain are 
summarized in Figure 34 as a function of dynamic pressure, 

LATERAL-DIRECTIONAL PERFORMANCE 

All preceding pitch axis time histories were obtained with the roll axis gain changer 
engaged. The variation in G^,p in response to pilot commands and turbulence exactly 
parallel the characteristics of the pitch gain The following summary will be 

limited to the performance of the lateral-directional CAS in conjunction with the gain 
changer. 

Yaw Axis Gain Changer Limitations 

Significant analysis effort was expended in attempting to apply the limit cycle to the 
yaw axis (rudder servo). The relatively low level of rudder effectiveness resulted in 
extremely high yaw rate gain as well as correspondingly high crossfeed gains. The 
yaw axis parameters — surface effectiveness, servo rate limit, servo loop gain, desired 
limit-cycle frequencies, and expected noise levels--were such that the loop could not 
be operated at critical gain levels without encountering excessive noise amplification 
resulting in severe servo rate saturation. This problem is not evident in pitch or roll 
due to higher combinations of surface effectiveness and servo rate. Consequently, the 
aileron servo was used as the basic input to the gain changer, and the yaw axis gains 
slaved to the aileron gain. 

Roll Axis Gain Changer Operation 

The gain changer in the roll axis is nearly identical to the pitch axis, and they probably 
could be made identical if any advantage would accrue. 

As in pitch, sensor noise will cause gain changer activity with the most significant gain 
reduction occurring for the assumed roll rate gyro noise. (During earlier yaw axis 
adaptive analysis, noise in all the yaw axis sensors and crossfeeds, coupled with high 
gains, would drive the adaptive gain to extremely low levels. ) Figure 35 is a summary 
of the adaptive gain level with yaw rate gyro, lateral accelerometer, roll rate gyro 
noise, and aileron servo hysteresis inputs as a function of flight condition. 


135 






0.020 

0.010 

0.005 

0,002 

0.001 

0.080 

0.060 

0.04 Q 

0.020 







Appendix E contains a set of time histories showing the response of the lateral-directional 
CAS with LCGC. The adaptive gains increase the damping of a step beta gust response 
from 0.15 (free aircraft) to about 0. 25 with little change in frequency. The turn 
coordination properties are satisfactory at all flight conditions. 



SECTION 11 

CONCEPT COMPARISON AND SELECTION 


Three factors were used to compare candidate adaptive concepts: 

• Performance, 

• Growth potential, and 

• Computer capacity requirements. 

On this basis the MLE concept was selected for further refinement and recommended 
for eventual flight test. Its strength lies primarily in growth potential. It also has a 
slight edge in performance relative to other concepts. 

PERFORMANCE 

The performance of each concept was evaluated against three "measures of goodness : 

. Accuracy at fixed flight conditions (including effects of pilot commands, 
turbulence, sensor noise, and actuator hysteresis), 

• Convergence properties, and 

. Tracking varying flight conditions. 

Tables 27 through 30 summarize the performance characteristics of each concept 
at four fixed flight conditions using a standard test sequence. For this test the initializa- 
tion was at the proper flight condition. The test sequence consisted of: 

1. Quiescent condition, 

2. Pilot doublet followed by four seconds, with no input, 

3. Thirty seconds of 99 percent turbulence (Dryden model), 

4. Pilot doublet with 99 percent turbulence (Dryden model), and 

5. Repeat of (1) through (4) with sensor noise. 

Flight conditions 1, 5, 8 and 10 are presented in Tables 27 through 30, respectively. 
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TABLE 27. PERFORMANCE CHARACTERISTICS AT FC1 
[20, 000 ft, M - 0.67, q = 14603 n/m 2 


1 ' ■ 

Tracker 

Limit cycle 

MLE 

Max. % errors in 

Excitation 

% error in g^ 

% of 

critical gain 

M 6 


V 

Kli 

Doublet 

5 

63 

4 

2 

6 


Turbulence 

t ; 

10 

75 

7 

2 

8 


Sensor noise 

20 

30 

4 

10 

15 

i 

■> Turbulence plus 
j sensor noise 

25 

30 

7 

2 

1 


Servo hysteresis 
with 0. 3 ft/s 
turbulence 

0 

74 

4 

2 

6 



a Angle-of-attack referenced to water-line. 


TABLE 28. PERFORMANCE CHARACTERISTICS AT FC5 
[20,000 ft, M = 0. 4, q = 5219 n/m 2 


i ' : 

Excitation 

Tracker 

M A 

% error in 

Limit cycle 
%of 

critical gain 

MLE 

Max. % errors in 

M 6 


B 

a a 
®L 

1 Doublet 

20 

75 

5 

H 

8 


S Turbulence 

20 

56 

5 

8 

2 


Sensor noise 

70 

20 

5 

B 

17 

8 

Turbulence plus 

20 

20 

5 

fl 



sensor noise 




■I 




Angle-of-attack referenced to water-line. 
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TABLE 29. PERFORMANCE CHARACTERISTICS AT FC8 
[40, 000 ft, M = 1. 2,q - 18913 n/m 2 



Tracker 

A 

% error in g^ 

Limit cycle 

MLE 

Max. % errors in 

Excitation 

% of 

critical gain 


M 

a 

H 

or T 
o»L 

Doublet 

30 

63 

0 

2 

■ 


Turbulence 

40 

81 

6 

6 

■ 


Sensor noise 

25 

43 

0 

2 

4 

4 

Turbulence plus 
sensor noise 

40 

43 

12 

2 

4 



TABLE 30. PERFORMANCE CHARACTERISTICS AT FC10 
[Sea level, M = 0.7, q = 34713 n/m 2 



Tracker 

A 

°Jo error in g 1 ^ 

Limit cycle 

MLE 

Max. % errors in 

Excitation 

% of 

critical gain 

M 5 

■ M 

a 

■ 

°u>L 

Doublet 

5 

48 

6 

9 

29 


Turbulence 

5 

69 

23 

9 

23 


Sensor noise. 

5 

> 48 

29 

5 

29 

3 

Turbulence plus 
sensor noise 

5 

48 

23 

9 

' \ 

23 



From these tables it is evident that sensor noise is mdre degrading to the tracker 
and limit-cycle concepts than to the MLE concept. MLE exhibits better accuracy 
than the model tracker. The limit-cycle design shows gain reductions from critical 
for pilot commands and turbulence as desired. • ' • , 
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The convergence characteristics of the three concepts are summarized in Table 31. 

The adaptation was determined for two types of inputs: 1) a sequence of pilot doublets 
(square wave C*-commands of 6.1 m/sec2 with a six- second period) with no turbulence, 
and 2) 99 percent turbulence with no pilot inputs. Two conditions were examined. For 
the first one, the aircraft was at FC1, but the algorithms were initialized at FC5. 

The MLE design converges faster than the model tracker. This test was not run for the 
limit cycle since it would imply initializing above critical gain which is not a realistic 
test. The second condition shows the aircraft at FC5 but initialized at a much higher 
dynamic pressure condition (FC10). Again, MLE has a faster response than the model 
tracker. During pilot commands, the limit-cycle design does not increase the loop 
gain to the proper value for this flight condition. During turbulence, the action of the 
set point does slowly return the gain to the appropriate value. 

TABLE 31. ADAPTATION CHARACTERISTICS WITH INITIALIZATION' 

AT DIFFERENT FLIGHT CONDITIONS 


Condition 

Input 

Time(Sec) to 80% Steady see 4 H 

Tracker) 

Limit cycle 

MLE 


Doublet 

4.0 

N. A. 

<1 

At FC1 

sequence 




initialized 

FC5 

99% 

turbulence 

4.0 

N. A. 

<i 


Doublet 

7.0 

- 

<1 

At FC5 

sequence 




initialized 

FC10 

99% 

turbulence 

7.0 

10 

<1 


Table 32 shows the tracking characteristics during an acceleration at maximum power 
from Flight Condition 5 (h = 6096 m, M = 6.4) to 7620 m and M = 1.1. Three 
different conditions (quiet, turbulence, and turbulence plus sensor noise) are presented. 
For the limit-cycle concept, it is desirable to maintain at least 50 percent of linear 
critical gain. For the two explicit concepts, it is desirable to accurately estimate 
surface effectiveness parameters. 
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TABLE 32. TRACKING CHARACTERISTICS FOR VARYING 
FLIGHT CONDITION 


Excitation 

Tracker a 
% error in gj ^ 

Limit cycle 
% of 

critical gain 

Peak MLE 

A 

% error in M c 
00 

Quiet case 

a 

25 

100 

a 31 

Turbulence 

20 

75 

28 

Turbulence 
plus sensor 
noise 

20 

36 

28 


a Test signal required. 


GROWTH POTENTIAL 

Growth potential is necessarily a qualitative notion. It is judged to be low for the 
limit-cycle system because the concept itself is restricted to a single (implicit) loop 
gain adjustment. 

In theory, the tracker provides greater growth potential because it can estimate several 
parameters as well as states. As discussed in Section 9, however, this potential is 
not realized in practice due to severe limitations imposed by sensor noise and disturbances. 
This leaves the MLE concept as the only approach with capability to achieve complex 
multidimensional gain adjustments. Its M gQ estimate can be used to adjust loop gain to 
compensate for changing surface effectiveness, while the M g estimate is useful for 
reduced-static-stability control laws. As discussed in Section 8, the JVT estimate can 
adjust an additional gain on a lagged pitch rate feedback to compensate for the changing 
stability properties of the airframe. This is especially important in transitions from 
subsonic to supersonic flight due to the shift in the center of pressure. 

COMPUTER REQUIREMENTS 

To estimate flight computer requirements, the adaptive subroutines were compiled 
and run on a CDC-6600. The resulting time and memory' requirements are listed in 
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Table 33. The data are given for FORTRAN generated code that is fairly efficient for 
the given algorithm (no indexing used). -A comparison of operation times for the CDC- 
6600 and the AP-101 flight computer is given in Table 34. This table suggests a speed 
ratio of 7:1 for typical multiply-add operations. Hence, the MLE control package which 
requires 1. 1 msec per control cycle on the CDC-6600 should consume 7. 7 msec per 
control cycle on the AP-101. The frame time for the F-8's digital system is 20 msec, 
so the control law plus adaptation will take well below half of the time available. However, 
this does not include time for redundancy management or mode switching which may be 
part of the total control repertoire. 

RECOMMENDED CONCEPT 

The concept comparison is summarized in Table 35. Overall, the MLE concept is 
judged to be the most promising and is recommended for refinement for flight test. 


TABLE 33. COMPUTER REQUIREMENTS BY SUBROUTINE 


Subroutine 

Memory 

CDC 6600 
Cycle time (msec) 

1* Pitch and lateral 

controller (no adaptive 
computation) 

570 (1072 octal) 

0.3 

2. Limit (up- and down- 
logic, both axes) 

546 (1042 octal) 

0.4 

3. Track 

272 (420 octal) 

0.2 

4. MLE 

2146 (4142 octal) 

0.8 
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TABLE 34. COMPUTER OPERATION SPEEDS 


Operation 3, 

CDC-6600 b 

IBM AP-101 C 

Memory access (i^sec) 

0.1 

0.9 

Add (\xsec) 

0.4 

2.4 

Multiply (nsec) 

0.7 

5.4 

Divide (p.sec) 

1.0 

10.0 - 10.5 


Floating point arithmetic, register to register. 

"L. 

D CD06600 CYBERNET applications group data. 
C IBM specification sheet data. 


TABLE 35. OVERALL COMPARISON OF CONCEPTS 


Quality 

Model tracker 

. 

■ 

High- gain 
limit cycle 

Maximum 

likelihood 

estimation 

Performance 

Acceptable to good 

Good 

Good to excellent 

Growth 

potential 

Low - limited to 
single variable 
explicit gain 
schedule 

Low - limited to 
single variable 
implicit adaptation 

High.,- multiple 
parameter gain 
adjustment 
possible 

Computer 

requirements 

Minimal 

Low 

----- _J 

Medium 


These requirements are relative to AP-101 capacity. 
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SECTION 12 


CONCLUSIONS AND RECOMMENDATIONS 


CONCLUSIONS 

Simulator quality designs for the F-8C DFBW application were successfully 
developed from three different adaptive concepts: maximum likelihood, model 
tracking, and high-gain limit-cycle control. Performance, growth potential, 
and computer requirements were used as criteria for selecting the most pro- 
mising of these candidates for further refinement. The maximum likelihood 
concept was selected primarily because it offers the greatest potential for 
identifying several aircraft parameters and, hence, for better control performance 
in future aircraft applications. In terms of identification and gain adjustment 
accuracy, the MLE concept proved slightly superior to the rest, but, as expect- 
ed, this increment has no significant effects on the control performance 
achievable with the F-8C aircraft. 

The MLE design is based on standard maximum likelihood estimation theory 
but uses a parallel channel implementation in order to avoid on-line data 
management and minimization iterations. It identifies a reduced set of three 
parameters from pitch axis data only. This structure was selected on the basis 
of separate theoretical identifiability studies which indicate that only limited 
identification potential exists for the F-8C under the ground rules imposed in 
Section 4. Identifiers which look for more parameters and/or use additional 
lateral-directional data cannot be expected to perform substantially better. Actual 
identification accuracies of the MLE concept are 10 to 20 percent for the critical 
M _ parameter. This is sufficient to schedule gains in both the pitch and 
lateral-directional control laws. 

The model- tracking design is based on modified Liapunov-stable model -following 
theory. Because this theory ignores process and measurement noise, it produces 
sensitive designs with large identification errors in applications which include 
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such disturbances. Hence, the F-8C model tracker relies on bandpass filters, thresholds, 
and limiters to minimize these, problems. Like MLE, the tracker was set up to identify 
three parameters. However, only the surface effectiveness estimate proved reasonable 
in the presence of gusts and sensor noise. 

The high-gain limit-cycle design is based on established flight tested control technology. 

It is applicable whenever the control structure has a single-loop gain as its variable 
parameter. Within this restriction, the design performed well. Its up-logic successfully 
alleviates the excessive gain reductions during pilot commands which were characteristic 
of previous flight systems, and there is no adverse coupling between simultaneous pitch 
and roll limit-cycle gain changers. 

RECOMMENDATIONS 

The maximum likelihood concept is recommended for flight test as part of the se cond 
phase of NASA's F-8C DFBW program. Several refinements of this design are also 
suggested. These include gust level estimation to fine-tune filter gains and automatic 
data length adjustment on the basis of aircraft maneuvers to reduce tracking errors. 

We further recommend development of a special purpose lateral identifier to extract 
improved velocity estimates during large rolling maneuvers. These estimates could 
be used to implement an inertially coordinated lateral CAS 2 without accurate air data. 
Advantages of inertial coordination are good dutch roll damping and good turn coordination 
at all angles -of-attack. These are not simultaneously achievable with the present reduced 
measurement lateral control law. 

Flight demonstrations of the MLE concept will mark an important step in the development 
of modern adaptive flight control technology. They will serve to verify quantitative 
predictions of identifier and control -loop performance and to establish qualitative 
properties of the overall control system. The latter include test signal requirements 
and their subjective effects on the pilot, performance during large maneuvers and 
configuration changes, the feasibility of quiescent operation without test signal excitation, 
and effects of modeling limitations associated with transonic flight, static and dynamic 
flexibility, real sensors, and actuators,, etc. Positive flight test results will contribute 
substantially to making adaptive control a viable alternative design approach for future 
aircraft applications. 
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APPENDIX A 


SCATTER PLOTS OP PITCH AXIS 
MODEL COEFFICIENTS 

x * Subsonic flight condition 

s * Supersonic flight condition 

Data does not include quasistatic 
aircraft flexibility* 
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Figure 56. Concluded. 
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Figure 57. MLE tracking test. Transition 1 
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Figure 57. Concluded. 
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Figure 58. MLE tracking test. Transition 1 in 99 percent turbulence. 
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Figure 58. Concluded. 
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Figure 59. MLE tracking test. Transition 1 in 99 percent 
turbulence plus sensor noise. 
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Figure 59. Concluded. 
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Figure 60. MLE tracking test. Transition 2. 
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Figure 60. Concluded. 
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Figure 61. MLE tracking test. Transition 2 in 99 percent turbulence. 



Figure 61 . Concluded. 


199 



Figure 62. MLE tracking test. Transition 2 in 99 percent 
turbulence plus sensor noise. 
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Figure 62. Concluded, 



Figure 63. MLE lateral maneuvers. 
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Figure 63. Concluded. 
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Figure 64. MLE roll over at FC1. 
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Figure 64. Concluded. 


205 





Figure 65. MLE a gusts 
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Figure 65. Concluded. 
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Figure 66. MLE gear-down transients. 
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Figure 66. Concluded. 
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Figure 67. MLE gear- up transient. 
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Figure 67, Concluded. 
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Figure 68. MLE wing-up transient. 
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Figure 68. Concluded. 
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Figure 69. Relaxed-static-stability flight, FC1. 
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Figure 70. MLE tracking test. FC1, no sensor noise ; 
half test signal. 
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Figure 70. Concluded. 
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APPENDIX C 


MODEL TRACKER TIME HISTORIES 


APPENDIX C 


MODEL TRACKER TIME HISTORIES 


The time histories taken on the Langley simulator consisted of 16 channels of data 
(two eight- channel recorders). The channels are: 


Recorder 1 
Pitch stick deflection 
Angle-of-attack (a) 

Pitch rate (q) 

Normal acceleration (a^) 
Elevator position (6^) 
Total velocity (V) 
Altitude (H) 

Flap position (6f) 


Recorder 2 
Dynamic pressure (q) 
q error 
Bank angle (cp) 


11 


12 


6 11 

G C* gain 

ET (output of nonlinear error 
function) 


Scales for each of these parameters are indicated in MKS and conventional units. 


ACCURACY AT FIXED FLIGHT CONDITIONS 


Time histories for standard pitch sequences consisting of periods of quiet, turbulence, 
and doublet commands are presented in Figures 71 through 75 for Flight Conditions 
1, 5, 8, 10 and 17. This sequence of runs was then repeated with gyro noise, 
accelerometer noise, and servo measurement noise simultaneously applied (Figures 76 
through 79). The sequence for Flight Condition 1 was repeated with servo flow hysteresis 
of 0.05 degree (Figure 80). 

CONVERGENCE RUNS 

Figure 81 shows the response when the model parameters are initialized at FC10 
but the aircraft is actually at FC5. Two separate inputs were applied; first turbulence 
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only and, secondly, a sequence of pilot doublets. Convergence rates are similar 
with time to reach 50 percent of steady state in two seconds, and 80 percent of steady 
state is reached in about four seconds. The gain response is slower due to the one- 

A 

second lag on the g^ parameter. 

Figure 82 shows a similar response when the model is initialized at FC5 and the 
aircraft is actually at FC1. Again, 80 percent of steady state is reached in about 
four seconds. 

PARAMETER TRACKING 

Two profiles corresponding to changing flight conditions were made. The first accelerated 
from FC5 (q = 105 psf) with full power (and afterburner). 

Figures 83 through 85 show the accelerating condition. For the quiet case, no adjust- 
ment in the model parameters occur and the errors grow as expected. Therefore, 
the quiet run was repeated using a random test signal bandlimited to have most of the 
energy around the short-period frequency. Use of a test signal held the tracking error 
to around 20 to 25 percent (Figure 83). For the case of turbulence, the tracking 
error in g^ is about 20 percent with a slightly large initial error buildup (Figure 84). 

Tracking during turbulence and sensor noise is shown in Figure 85. 

The second profile consisted of a deceleration from Flight Condition 10. Thrust is 
reduced to idle and the speed brake extended. Airspeed gradually decreases and altitude 
remains essentially constant. The quiet condition (Figure 86) again required a test 
signal. Peak errors are on the order of 50 percent. In tracking during turbulence 
(Figure 87), the error in g^ is less than 20 percent. For the turbulence and sensor 
noise condition, the peak errors are 20 percent (Figure 88). 

LATERAL DISTURBANCE 

Of interest is the effect of a roll maneuver on the pitch axis tracker. Figure 89 shows 
the time history of a snap 360° roll at Flight Condition 1. This is a severe test for the 
tracker since the sensor bandpass filters cannot eliminate the variation in the normal 
accelerometer reading. The peak error developed is 60 percent. Slower roll rate will 
have less effect. The figure shows that any activity following the roll will reduce the error. 
In the example, turbulence is applied following the roll and the error is reduced. 
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Figure 71. Model tracker at FC1. 
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Figure 71. Concluded. 
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Figure 72 


Model tracker at FC5 



Figure 72. Concluded. 
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Figure 73. Model tracker at FC8. 
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Figure 73. Concluded. 
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Figure 74. Model tracker at FC10. 
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Figure 74. Concluded. 
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Figure 75. Model tracker at FC17 
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Figure 75, Concluded. 
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Figure 76. Model tracker at FC1 (sensor noise). 
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Figure 7.6. Concluded. 
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Figure 77. 


Model tracker at FC5 (sensor noise). 
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Figure 77. Concluded. 
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Figure 78, 


Model tracker at FC8 (sensor noise) 



Figure 78. Concluded, 
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Figure 79 


Model tracker at FC10 (sensor noise) 



Figure 79. Concluded. 
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Figure 80. Model tracker at FC1 (with hysteresis). 
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Figure 80. Concluded. 
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Figure 81. Convergence from FC10 to FC5. 
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Figure 81. Concluded, 
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Figure 82* Convergence from FC5 to FC1. 
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Figure 82. Concluded. 
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Figure 83. Model tracker during acceleration 
(test signal). 
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Figure 83. Concluded. 
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Figure 84, Model tracker during acceleration (gusts), 



Figure 84. Concluded. 
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Figure 85. Model tracker during acceleration 
(gusts plus sensor noise). 
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Figure 85. Concluded. 
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Figure 86. Model tracker in deceleration (test signal). 
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Figure 86. Concluded. 
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Figure 87. Model tracker in deceleration (gusts). 
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Figure 87. Concluded. 
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Figure 88. Model tracker in deceleration 
(gusts plus sensor noise). 
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Figure 88, Concluded. 
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Figure 89, Model tracker response to snap roll (FC1). 
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Figure 89. Concluded. 
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APPENDIX D 


LIMIT-CYCLE GAIN- CHANGER 
PITCH AXIS TIME HISTORIES 


APPENDIX D 


LIMIT- CYCLE GAIN- CHANGER 
PITCH AXIS TIME HISTORIES 


Time histories taken on the Langley simulator consisted of 16 channels of data (two 
eight-channel recorders )• The channels are: 


Recorder 1 
Pitch stick deflection 
Angle-of-attack (a) 

Pitch rate (q) 

Normal acceleration (a 3 > 
Elevator position (6^) 
Total velocity (V) 
Altitude (H) 

Flap position (6^) 


Recorder 2 
Dynamic pressure (q) 

G C* g ain 

Elevator servo position (6 ) 

e s 

Aileron servo position (6 &g ) 
G LAT gain 

XX 

nk angle (cp) 

Roll rate (r) 

Sidelsip (0) 


Scales for each trace are indicated in MKS and conventional units. 


ACCURACY AT FIXED FLIGHT CONDITIONS 


Standard pitch sequences consisting of periods of quiet, turbulence, and doublet 
commands are presented in Figure 90 through 94 for Flight Conditions 1, 5, 8, 10 
and 17 (power approach). This sequence of runs was then repeated with gyro noise, 
accelerometer noise, and servo measurement noise simultaneously applied (Figure 95 
through 98). 

Finally, the sequence at Flight Condition 1 was repeated with servo hysteresis of 0.05 
degrees. (Refer to Figure 99.) 
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CONVERGENCE PROPERTIES 


The concept of "converging" to a gain value does not apply to the limit-cycle concept 
in the same way it does to ttrfe explicit identifiers. In normal operation, the gain 
changer will not be at a gain higher than critical. If the system starts at a value below 
critical, the gain increase is largely determined by the integral gain. Figures 100 
and 101 show the responses at FC5 with the system initialized at FC10. In turbulence 
it takes on the order of ten seconds to reach "steady state" (which is six db below 
critical). (Refer to Figure 100). For pilot commands, the gain is reduced slightly. 

ADAPTATION 

In the case of increasing and decreasing dynamic pressure, the gain changer tracked 
the critical gain (quiet case) or tracked the "turbulence" gain value. Figure 102 
shows a full power acceleration starting from FC 5. Note the constant limit-cycle 
amplitude. Figures 103 and 104 show the same run with turbulence and then with 
turbulence and sensor noise. 

Figures 105 through 107 show the responses starting from FC 10 and decelerating by 
reducing thrust to idle. Again the limit-cycle amplitude remains constant (Figure 105). 
Adaptation during turbulence is shown in Figure 106. The same adaptation is repeated 
in Figure 107 with turbulence and sensor noise applied. 
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Figure 90. Limit-cycle design at FC1. 
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Figure 90. Concluded. 
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Figure 91. Limit-cycle design at FC5. 


264 



Figure 91. Concluded. 
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Figure 92. Limit-cycle design at FC8. 
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Figure 92. Concluded, 



Figure 93. Limit-cycle design at FC10. 
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Figure 93. Concluded. 
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Figure 94. Limit-cycle design at FC17. 
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Figure 94 
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Figure 95. Limit-cycle design at FC1 
(sensor noise). 
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Figure 95. Concluded. 
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Figure 96. Limit cycle design at FC5 (sensor noise). 


274 



Figure 96. Concluded. 
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Figure 97 


Limit-cycle design at FC8 (sensor noise) 
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Figure 98. Limit-cycle design at FC10 (sensor noise) 



Figure 98. Concluded. 
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Figure 99. Limit-cycle design at FC1 

(sensor noise plus hysteresis). 
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Figure 99 • Concluded. 
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Figure 100, Limit-cycle design at FC5, IC FC10, gusts. 
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Figure 100. Concluded. 
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Figure 101. Limit-cycle at FC5, IC FC10, pilot commands. 
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Figure 102. Limit-cy6le design under acceleration. 
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Figure 102. Concluded. 
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Figure 103. 


Limit- cycle design under acceleration (gusts). 
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Figure 103. Concluded. 
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Figure 104. Limit-cycle design under acceleration 
(gusts plus sensor noise). 
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Figure' 104. Concluded. 


290 


<1 •> 
n/in 
(psf) 


G A 


6 . 

o 

s 


rad 



G 


I .AT 


V 

rad 


P 

rad /see 
(deg/sec) 


& 

rad 

(deg) 



Figure 105. Limit-cycle design, deceleration from FC10. 
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Figure 105. Concluded. 
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Figure 106. Limit-cycle design deceleration from FC10 (gusts). 
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Figure 106. Concluded. 
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Figure 107. Limit-cycle design deceleration from FC10 
(gusts plus sensor noise). 
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Figure 107. Concluded. 
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APPENDIX E 


LIMIT- CYCLE GAIN-CHANGER 
LATERAL AXIS TIME HISTORIES 


The series of analog traces (Figures 108 through 118) shows the response of the 
lateral-directional CAS to beta gusts and roll rate commands. The adaptive gains 
increase damping of the aircraft beta gust response from about 0. 15 to about 0. 25, 
with little or no change in frequency. In contrast, the nominal reduced CAS of 
Reference 1 tended to normalize the response frequency to about 0. 4 Hz with similar 
damping. The tendency of the adaptive gains to reduce to the six-db level for command 
inputs is evident. 

For further comparison, the gain was fixed at six db lower than critical, and the 
response is nearly identical to the adaptive response. The responses to lower gains 
are also shown since under turbulence and roll rate gyro noise, lower gain levels are 
realized with the gain changer. This data was also used to assess the applicability 
of using a gain scheduled on estimated surface effectiveness (or dynamic pressure) 
for use with the two other adaptive concepts. The performance degradation of 
operating the roll axis at a lower gain is minimal. The yaw axis gain does not need 
to be reduced in a scheduled system since it has sufficient gain margin. 

Figure 119 shows roll rate commands and beta gust responses from the Langley 
F-8C simulation. Flight conditions 1, 5, 8, 10 and 17 are presented. They compare 
well to the analog responses for the corresponding flight condition. 
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Figure 108. Flight Condition 1. 
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Figure 108. Concluded. 
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Figure 109. Flight Condition 5. 
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Figure 109, Concluded. 
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Figure 110. Flight Condition 6. 
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Figure 110. Concluded. 
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Figure 111. Flight Condition 7. 
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Figure 111. Concluded. 
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Figure 112. Flight Condition 8. 
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Figure 112. Concluded. 
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Figure 113. Flight Condition 9. 
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Figure 113. Concluded. 
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Figure 114. Flight Condition 10. 
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Figure 114. Concluded. 
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Figure 115. Flight Condition 11. 
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Figure 115. Concluded. 
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Figure 116. Flight Condition 12. 
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Figure 116. Concluded. 
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Figure 117. Flight Condition 16. 
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Figure 117. Concluded. 
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Figure 118, Flight Condition 17. 
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Figure 119. Lateral Responses from LRC F-8C Simulation 
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